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FACTORS CONCERNED IN THE PHENOMENA OF OXIDA- 
TION AND REAERATION? 


INTRODUCTION 


The studies presented in this report were made by the United 
States Public Health Service during the years 1914 and 1915 as part 
of a comprehensive survey of the pollution and natural purification 
of the Ohio River, conducted under the supervision of Surg. W. H. 
Frost. The scope and purposes of this survey as a whole are out- 
lined in a previous publication,’ to which reference is made for a 
detailed description of the Ohio River, summaries of its sources of 
- pollution, and measurements of discharge and velocity, also for 
presentation and discussion of the results of bacteriological exami- 
nations and chemical analyses other than those dealing with observa- 
tions on dissolved oxygen. These latter observations, which are 
presented in the previous publication only in the form of a basic 
summary * without discussion, are regarded as being sufficiently dis- 
tinctive in their significance and in the character of analysis required 
to justify their separate treatment in this supplementary report. 

In the case of a stream like the Ohio River, important not only for 
its size and navigability but also because it is the sole available source 
of water supply for a large and growing population, permissible 
limits of its pollution are given in bacteriological rather than chemi- 
eal terms; that is, they are fixed by the permissible load of bacterial 
pollution which may be placed upon the river with reference to water 
supplies rather than by considerations involving the exhaustion of 
the oxygen supply and physical nuisance. Nevertheless, the possi- 
bilities existing for at least a partial depletion of the reserve oxygen 
supply of the river, under conditions of pollution already approach- 
ing a critical stage from a bacterial standpoint in certain zones, could 
hardly be neglected in any consideration of the general problem. 
An exceptional opportunity was presented, moreover, for a study 
of the oxygen status of the river, and especially of its purification 
capacity from an oxygen standpoint, in view of the availability of 
extensive collateral data relative to river temperature, discharge, 


1 Manuscript submitted June, 1924. 

2A Study of the Pollution and Natural Purification of the Ohio River. II. Report on 
Surveys and Laboratory Studies, Public Health Bulletin No. 143, Washington, 1924. 

3 Op. cit., Table No. 51, pp. 124-128. 
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and times of flow between the various sampling stations, collected 
in connection with the investigation. 

The general viewpoint of the studies here described was influenced 
to a large extent by the newer conception of stream purification 
which has resulted from the marked advances that have taken place 
during recent years in the application of bacteriology and of physical 
and biological chemistry to sanitary science. First: in importance 
have been discoveries relative to the generality of laws governing 
the death rates of bacteria, which have emphasized the progressive 
character of the complex biochemical reactions concerned in stream 
purification, and hence the controlling influence exerted by the time 
factor in such phenomena. Of almost equal significance has been 
the evolution of a newer biochemistry of sewage and sewage-pol- 
luted waters, wherein the older sewage chemistry, dealing with 
nitrogen in its various forms, has been largely replaced, in prob- 
lems involving the stability of organic matters of sewage origin, by 
biochemical methods of study permitting a direct measurement of 
the oxidation reactions more directly related to organic stabiliza- 
tion processes. As examples of these methods may be cited the 
“relative stability ” and “ biochemical oxygen demand ” tests, which 
are familiar to everyone who has followed the literature of sewage 
during recent years. Finally, the great value of modern physical 
chemistry as an aid in interpreting and applying to stream condi- 
tions the results of biochemical methods of study should be noted. 

Previous to the foregoing developments, studies of the self-puri- 
fication of streams had necessarily been, in the main, empirical; that 
is, they had comprised. the determination and recording of actual 
conditions measured in analytical terms, without any attempt to 
formulate results in terms of general principles. Such a procedure 
is valuable as a matter of historical record relative to a given stream 
_ or local condition, but unfortunately it fails to give data having pos- 
sibilities of more general application. In some cases generalizations 
such as have been attempted from data of this character have re- 
sulted in serious misconceptions of the relative importance of dilu- 
tion as a factor in the oxidation of waste matters in streams, as 
compared with that of reaeration, which is often of far greater 
significance. A simple example will illustrate this point. 

Accepted standards for a safe dilution ratio, based on certain 
stream conditions, especially in Massachusetts, are from a minimum 
of 3.5 to a maximum of 6.0 second-feet of normal unpolluted stream 
water per thousand of population contributing sewage. 


The sewage of the District of Columbia has in summer a total bio-_ 


chemical oxygen demand equivalent to 112 grams per capita daily, 
or about 300 parts per million when corrected to a normal sewage 
strength basis of 100 gallons per capita flow, which is equivalent to 
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roughly 0.15 second-feet per thousand of population. These’ oxygen 
demand values, which agree very closely with the figure, 100 grams 
per capita daily, given by Pearse‘ as a result of studies of Chicago 
sewage, may be taken as being fairly Hol gas BML for noe 
domestic sewage. 

During the summer period, June 1 to October 15, 1914, the mean 
discharge of the Ohio River at a point immediately bate Cideilati 
represented a flow of 5.1 second-feet per thousand of urban popula- 
tion in the watershed above this point, at which point the average 
dissolved oxygen content observed was 5.7 parts per million, or 68 
per cent saturation. Assuming the normal summer oxygen content 
of a theoretically unpolluted Ohio River to be at the mean summer 
saturation value, 8.2 parts per million, the draft imposed upon the 
normal dissolved oxygen content of the river to satisfy the oxygen 
demand of the urban sewage alone would be represented, in this 
case, by the difference between 8.2 and 5.7 parts, or 2.5 parts per 
million. On this basis, the required dilution, in the absence of re- 
aeration, would be 300: 2.5, or 120: 1, amounting to a stream flow of 
18.5 second-feet per thousand of population. Since the flow was 
actually but 5.1 second-feet per thousand, it would appear. that 
dilution alone was responsible for about one-quarter, and reaeration 
for nearly three-quarters, of the total purification of the river taking 
place up to the point in question. With a permissible oxygen content 
of the Ohio River at this point lower than the approximately 70 per 
cent saturation figure observed during the summer of 1914, the stream 
flow actually required would be reduced by even more than a propor- 
tionate amount, as the relative quantities of dissolved oxygen sup- 
plied by reaeration would be greater with lower saturation values in 
the river. | 

The foregoing example not only illustrates the great importance 
of reaeration in the purification of running streams, but also em- 
phasizes the fallacy of depending solely upon mere dilution as a 
measure of the oxidation of wastes discharged into flowing bodies 
of water. An attempt to apply dilution ratios, derived from rivers 
having great reaeration capacities, to sluggish streams like the 
Chicago drainage canal, would unquestionably lead to serious error, 
as has been borne out by experience with a number of sluggish canals, 
which have given rise to offensive conditions, though not excessively 
polluted from a standpoint of ordinary dilution criteria. 

The development of a better understanding of the relations be- 
tween dilution and reaeration as factors in stream purification, and 
in fact our modern conception of the overwhelming importance of 
dissolved oxygen in determining the power of natural bodies of 
water to digest and oxidize organic polluting matters, owes much to 


4Trans. Am. Soc. C. E., vol. 85, 1922, p. 451. 
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the. pioneer work in sewage bio-chemistry carried out by the Frank 
lands, McGowan, Letts, Dibdin, Adeney and their coworkers ir 
England, and to adaptations and modifications of their methods o! 
study by various workers in the United States, notably Black anc 
Phelps, Hoover, Lederer, and Theriault. While no extended bibliog: 
raphy of the studies above noted will be presented in connectior 
with the present paper, a number of references covering certair 
specific points will be given in the later text. . . 

With the aid of the newer bio-chemistry of sewage and of the 
quantitative view of the oxygen relations in polluted streams whick 
it permits, it has been possible to formulate and test empirically ¢ 
general theory of stream purification from an oxygen standpoint 
using for this purpose the data obtained in connection with the Ohic 
River studies. Before proceeding with a presentation and analysis 
of these data, it will clarify the subsequent discussion to outline this 
underlying theory. 


THEORETICAL DISCUSSION 


_ The capacity of a stream to receive and oxidize sewage depend: 

upon its, oxygen resources. The condition of a polluted stream a 
any time is the result of a balance between these resources and th 
demand made upon them by the organic polluting matter carried by 
the stream. This demand, being the result of a slow bio-chemica 
reaction, is, in the absence of new pollution, a progressively decreas 
ing one, and as the resources of the stream are composed in part of : 
continuous influx of oxygen from the atmosphere, the state of bal 
ance which determines the momentary condition of the stream 1 
constantly changing. There are, therefore, two primary phases 0. 
the problem—namely, the actual, momentary condition, and th 
direction and extent of the existing changes, which indicate th 
future condition. Fresh sewage, for example, may contain some dis 
solved oxygen, and, measured upon the oxygen scale of nuisance 
may be in the same momentary condition as a stream which has abou 
completed the work of oxidizing organic pollution and contains th 
same amount of residual dissolved oxygen. The direction of change 
however, is entirely different and determines the distinction betwee 
the two cases. The oxygen resources are comparable to the assets o 
a balance sheet; the oxygen demand to the liabilities. The condition 
of a strong sewage containing oxygen is comparable, in financia 
terms, to one of momentary solvency, with available cash, but witl 
excessive obligations maturing on the morrow. A comprehensiv 
theory of self-purification must therefore deal with the oxyger 
demand as well as with the oxygen resources, and must consider th 
relation of the various factors of time, temperature and other phys: 
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cal conditions to the rates of change of these two fundamental ane 
tities. 
THE OXYGEN DEMAND UPON A STREAM 


Changes in the dissolved oxygen content of a stream are inti- 
mately associated with biochemical changes. They are brought 
about primarily by the oxidation of organic matter discharged into 
streams as soil wash and as wastes. In the presence of a supply of 
oxygen, together with certain oxidizing bacteria and oxidizable or- 
ganic matter, progressive oxidation and stabilizing of the organic 
matter will take place. 

It has been shown? that, under experimental conditions approxi- 
mating those prevailing in a stream containing reserve dissolved 
oxygen, this reaction is an orderly and consistent one, prececame at 
a measurable rate and according to the following eens law :? 

The rate of biochemical oxidation of organic matter 1s propor- 
tional to the remaining concentration of unoxidized substance, meas- 
ured in terms of oxidizability. 

This law, which happens to be similar to that which defines the 
course of a monomolecular reaction,’ may be stated in its differential 


form thus: 


dL K 
a 
a 


which may be integrated to the form— 


logy =Kt 


L’ being the initial and L the final oxidizability, or oxygen demand 
of the organic substance, in terms of oxygen; t being the elapsed 
time and K a constant coefficient, defining the rate at which the re- 
action proceeds. The value of K depends upon the character of 
organic matter and upon the temperature. 

It will be assumed for the present that this relation holds dd 
under actual stream conditions. Evidence supporting this assump- 
tion will be presented later in the text (p. 40). 3 

Defining the oxygen demand as the total remaining oxidizability 
of the substance present at any time, the law states that in equal 


1 Phelps, Harle B., Biochemistry of Sewage, VIII, Int. Cong. Appl. Chem., XXVI, 251. 

2This statement should be qualified to the extent of noting that little definite knowl- 
edge exists as to whether the law stated holds for periods of time longer than about 20 
days. Experimental data bearing on this point are, in fact, somewhat meager for periods 
longer than 10 days, though for shorter periods the most reliable evidence has been 
confirmatory. 

’The similarity of this law to that of monomolecular chemical reactions is ree 
due to the biochemical nature of the phenomenon. If the reaction were strictly chemical 
it would theoretically be more likely to follow the law of bimolecular reactions, as the 
two reacting substances, oxygen and oxidizable organic matter, are both present in 
limited amounts, 
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periods of time an equal proportion of the remaining oxygen de- 
mand will be satisfied. That is, if 20 per cent of the initial oxygen 
demand be satisfied in the first 24 hours, 20 per cent of the remain- 
ing demand will be satisfied in the second 24 hours, and so on. 
(See fig. 1. for graphical illustration.) Simce the oxygen demand, 
as actually determined in the laboratory, is given in terms of 


Fig. No. i 
CHART SHOWING RELATION BETWEEN OXYGEN PEMAND 
AND TIME AT 20 ° CENTIGRADE 
Demand in ony time = 24 Hpur Demond x Be 
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with K,.=0.1, ond of temperature of BO °C 
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dissolved oxygen, the rate of satisfaction of the demand, which is 
denoted by the term (- ony is equivalent to the rate of oxygen de- 


pletion. If the dissolved oxygen content be expressed in terms of 


| Meee eo ; 
saturation deficit (D,), its rate of depletion is (Se) and the follow- 


ing relation therefore holds: 


The nature of the substances found in city sewage and the num. 
ber and kinds of bacteria present have been found to be sufficiently 


losmnd 
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constant under various conditions to give a fairly constant value of K 
for this reaction at a constant temperature. Its value at 20° C., 
the time being expressed in days, was determined by Phelps by the 
methylene blue reaction for the sewage of Boston,* and later by 
direct determination of the oxygen demand of the sewage of Wash- 
ington; and also computed from results obtained by Lederer’ at 
Chicago in a comparison of the dilution and saltpeter incubation 
methods, using methylene blue. From all these data, values of K 
in the above equation of approximately 0.1 were obtained or de- 
rived. Later work by Lederer and collaborators® has indicated 
some variation of value of K for different sewages, though these 
differences are not striking. Values of the constant for different 
organic industrial wastes appear to vary more, doubtless because 
of the widely varying character of the organic matter present in 
such wastes, though the results of recent studies by Theriault’ indi- 
cate that when tests are carefully controlled values of the constant 
similar to those for sewage are obtained for a wide class of industrial 
wastes. Values for polluted river water should not vary greatly 
from those for average sewage, since a major part of the oxidizable 
organic matter present in such water has its origin in sewage.® 

For any particular sewage or similar waste the value of the 
velocity constant K is a function of the temperature. The tempera- 
ture relation has been found empirically to be defined by the 


formula: 
ee) 


an 
in which (T’) and (T) are the two temperatures, (K’) and. (K) 
the corresponding values of the velocity constant of the reaction, 
and © the thermal coefficient, a constant for the reaction, which may 
be determined experimentally. — 7 

The value of the thermal coefficient © has been determined experi- 


mentally at Cincinnati in connection with the present studies and 
also by Phelps in connection with studies at Boston and at the 


4 Phelps, Earle B., U. 8. Geol. Survey, W. S. Paper No. 229, Washington, 1909. 

6 Lederer, A., Jour. Inf, Dis., Vol. XIV, 1914, 482. 

® Lederer, A., Am. Jour. Publ. Health, Vol. V, 1915, p. 354. 

7 Public Health Reports, Reprint No. 594 (1920), U. S. Public Health Service, Wash- 
meton, D.-C, 

§ Later studies of the value of K, made on waters of the Chicago drainage canal, the 
upper Illinois River, and the Ohio River at Cincinnati, have given results closely con- 
firming this statement. The studies on Ohio River water, which were in progress at 
the time of final revision of this text, have involved probably the most elaborate and 
earefully controlled series of tests thus far made of the form and constants of the oxi- 
dation curve, employing three separate temperatures of incubation, 8° C., 20° C., and 
30° C. For periods of time up to about 20 days these curves have been found to follow 
closely the law stated on page 5, and their (K) values at 20° C. to agree within from 
5 to 10 per cent with the value, 0.1, employed in connection with the present text. 
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Hygienic Laboratory at Washington. These values are given below, 
as follows: | : 


1. Cincinnati (Mill Creek water) : Value of © 
Temperature 'ranze, 10° Coto20® Cures seo get ea ee 1, 0524 
Temperature range, 20° -C: 101876 wes 44432 oat Bi a 1, 0504 

2. Boston (sewage, methylene blue method) : 

Temperature range, 20° C. to 87.5° C____--- Ye ee petigaseiens 1, 0415 

3. Hygienic Laboratory (nitrate method) : _ 
Temperature range)-15° C. to' 24.5". C2252 Sars PILROUE IAS oe 1. 0441 
Température range, 24)5° "Coto3% * Ogg 12 pt ssipes_sugyede * 1/01 

Average, excluding value marked (*)~~----__--__--------~---------- 1. 047 


For the purpose of the present discussion we will adopt the mean 
value of 1.047, obtained from the data given in this table, after 
eliminating the value of 1.01, which 1s obviously out of harmony 
with the others. The variation in the value of (K) at different tem- 
peratures, as defined by a value of © equal to 1.047 is illustrated 
in Figure No. 2, by the line showing the relation of K values at 
various temperatures to its value at 20° ©. Conversion of the K 
values from one temperature to another may be readily made by 
means of a similar chart. 

As stated on page 7 (footnote 8), further studies of the tem- 
perature relations of (K) for Ohio River water at Cincinnati, 
Ohio, were in progress at the Cincinnati laboratory of the Public 
Health Service at the time of final revision of the present text. From 
the results of these studies available up to that time, it appeared 
that values of the thermal correction factor © would be obtained 
agreeing closely with the mean value, 1.047, as above derived, assum- 
ing the experimental values of (K) at the different incubation tem- 
peratures to be computed in the same manner as in the present case. 
The method of computation thus followed has consisted of calecu- 
lating, for the oxidation curve obtained at each respective incubation 
temperature, the percentage of the 90-day oxygen demand at that 
temperature remaining at successive intervals of time up to 20 days. 
These values, when plotted as ordinates against incubation times, 
using semilogarithmic paper, define a straight line having a slope 
equal to (K). This procedure carries an assumption that, for all 
practical purposes, the oxygen demand satisfied in 20 days is equiv- 
alent to the total available oxygen demand (with a value of K at 
20° C. equal to 0.1, it is actually 99 per cent of the total). This 
assumption is substantially correct for temperatures of 20° C. and 
upwards, but becomes less so as the temperature is diminished below 
20° C., since at the lower temperatures the time required for the satis- 
faction of a large proportion of the total oxygen demand is longer 
than at 20 degrees. When the studies described in the present text 
were made, no basis existed for extrapolating oxidation curves 


9 


beyond a time of 20 days at any temperature; hence, the coefficient 
(K), whether uncorrected at 20° C. or corrected to its equivalents at 
other temperatures, can hardly be regarded as being well defined 
for any period of time exceeding 20 days. This limiting time, how- 
ever, usually is long enough for ordinary purposes of calculation. 
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THE OXYGEN RESOURCES OF A STREAM 


Unpolluted water tends to hold in solution the maximum amount 
of oxygen which it is capable of containing at the existing tempera- 
ture and partial oxygen pressure of the atmosphere. This is the 
so-called saturation value ® and ranges at normal sea level barometric 
pressure from approximately 14 mg. per liter at just above freezing 
to about 7.6 mg. per liter at 30°C. 

In polluted streams the draft imposed upon the dissolved oxygerm 
supply by the progressive satisfaction of the oxygen demand reduces: 


® The actual saturation value depends upon the partial solution pressure of oxygen im 
the overlying atmosphere, and is directly proportional to the percentage of oxygen in the 
overlying gas. For normal atmospheric air, the solution pressure, and therefore the 
oxygen saturation value, is approximately one-fifth its value with an overlying atmos-~ 
phere of pure oxygen. 


sl ae 

the oxygen content below its saturation value; but as soon as thi 
depression occurs absorption of oxygen from the atmosphere fol 
lows. The sources of available oxygen in a stream are therefor 
mainly two—the oxygen initially dissolved in the water, with it 
maximum possible equilibrium value at the saturation point, and th 
oxygen absorbed from the atmosphere by partially deoxygenate 
water in the process of physical reaeration.”° Reaeration is a facto 
of so great importance in the biochemical oxidation of a stream tha 
it merits discussion in some detail. | | 

In a quiescent, unsaturated body of water the surface film, i 
immediate contact with the atmosphere, rapidly becomes saturate 
with oxygen, and can dissolve no more of the gas until that whicel 
has already passed into solution at the surface has been carried t 
the deeper strata by diffusion or “streaming,” or by both of thes 
agencies comkined. A more or less well defined concentratio. 
gradient is thus formed, with approximate saturation at the surfac 
and decreasing concentrations at increasing depths, the draft on th 
surface film being balanced by replenishment of the supply throug! 
solution from the atmosphere. Thus the rate at which a give 
small volume of water becomes reoxygenated is a function of it 
saturation deficit and its distance from the surface. 

In a flowing stream conditions of quiescence are, of course, modi 
fied greatly by the constant vertical mixing or overturning whic! 
the water receives as it flows along the channel. An approach t 
quiescence may be found in occasional pools during periods of slug 
gish flow, but even under these circumstances slight currents an 
eddies, due to movements of the water and to slight differences 
temperature, bring about a more or less continuous mixing of th 
stream throughout its cross section. The concentration gradien 
characteristic of quiescent waters is therefore subject to constan 
disturbance, but since a comparatively small proportion of the wate 
reaches the exact surface layer, the forces distributing the dissolve 
oxygen throughout the body of the stream are those of mechanica 
mixing and diffusion combined. Viewed as a diffusion phenomenor 
this distributing action may be conceived as being one in which th 
effect of mixing is simply to greatly decrease the effective dept. 
through which diffusion acts, the amount of contraction in the effec 
tive depth being a function of the turbulence of the stream and de 





ee ES ee eee UE Ore! Sayer. 

10 A third possible source of available oxygen is the activity of certain aquatic plant 
notably the green and blue-green. algae, which occasionally find in sluggish pools a favo 
able environment for their propagation. This, however, is a special case which is ne 
ordinarily an important factor in the reaeration of flowing streams. 

The term “ streaming” has been applied by Adeney to the downward movement ¢ 
currents of water from the surface to the lower strata, due largely, as subsequently note 
by Phelps (Trans. Am. Soc. C. E., vol. 76, p. 1626), to slight differences in density resul 


ing from surface evaporation. ay 
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pending upon the various physical characteristics of channel and 
tlow effecting turbulence. 

The conception that diffusion is the active agent in the distribution 
of dissolved oxygen throughout a turbulent stream, in the presence 
of what might appear to be the overwhelmingly superior forces of 
mechanical mixing, is a rather novel one and perhaps not readily 
obvious. Mechanical mixing of the water constantly breaks up the 
saturated surface film into small masses of water which are brought 
into contact with less saturated masses with an intimacy and fre- 
quency depending upon the degree of agitation. But the oxygen 
cannot pass from one zone to another, no matter how thin the zone or 
intimate the mixture of particles, except by diffusion. If this force 
were absent, the effect of mixing, whether mechanically or by 
“streaming,” would be similar to that of combining a measure of 
peas and beans, the aggregation being merely one of saturated and 
unsaturated particles of water. The laws. of diffusion therefore 
have an important bearing upon the phenomenon of reaeration, even 
in running streams. 7 

The process of diffusion had been studied by many of the earlier 
chemists, but the first attempt to formulate a definite principle was 
made by Fick,!? who enunciated what has since been known as Fick’s 
law of hydrodiffusion, as follows: 

A dissolved substance tends to diffuse through the solvent from a region of 
greater to one of less concentration, and the amount of substance which passes 
through a unit cross section in unit time is proportional to the difference in 
concentration on either side of this section. 

This law in its general form is capable of expression only as a 
differential equation, which requires for its solution the establish- 
ment of certain restricting or boundary conditions. ' 

A mathematical analysis of the law of Fick was made by Black 
and Phelps** for the following special case. A column of water of 
stated depth initially contains dissolved oxygen at a stated concen- 
tration, less than that of saturation, and uniformly distributed 
throughout the length of the column. At the initial time this col- 
umn is exposed at its surface to the atmosphere the process of down- 
ward diffusion is inaugurated and a concentration gradient estab- 
lished from surface to bottom. Required, the average concentration 
of oxygen throughout the column at the end of a stated time. 

The final integrated equation which was obtained is— 


B og e7 8k e@7 25K e749 
D=100-(1~495 81.08 )(e 49-455 +e ---) 


12 Poge, Ann., 94, 1855; p. 59. Nernst, Theoretische Chemie, 2d German ed., p. 157. _ 
138 Black, Col. W. M., and Phelps, H. B., The Discharge of Sewage into New York Har- 
bor. Report made to the Board of Hstimate and Apportionment, New York City, 1911. 
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in which— 
D=the required average concentration of dieblved oxygen ex- 
pressed in per cent of saturation, aiter time, t. 
B=the uniformly distributed initial iicaiesaie expressed in 
the same units. 


e=2;718. 
Tester 
k= rial Ee 


t—the elapsed time since concentration (B) existed. 
L=the linear depth of the column. 
athe diffusion coefficient, a constant for pure water, and a given 
gas (oxygen) at constant temperature. 


This equation applies primarily to quiescent water, but it has 
already been indicated that the effect of mixing is virtually to decrease 
the effective depth through which diffusion acts, so that the applica- 
tion to turbulent waters is permissible, the depth term being then 
indeterminate. | 

For the purposes of the present discussion it is desired to determine 
- the rate of solution of oxygen under any stated conditions. This rate 
is derivable from the equation of condition by differentiation with 
respect to a variable time. The actual operation need not be per- 
formed, for it is obvious that the result will be a complex expression 


containing the factor (1 ey 00) the initial saturation deficit, together 


with a complex time function. The rate of solution, therefore, is 
proportional, among other things, to the initial saturation deficit, and 
the proportionality factor is itself a time function—i. e., it varies 
with the time. At zero time—i. e., at the start— the rate is strictly 
proportional to the deficit. 

In any stream having sufficient turbulence to make reaeration a 
real factor in self-purification, there is sufficient mixing to bring 
about a fairly uniform condition throughout the vertical cross- 
section, so that at no time is there material vertical stratification. 
Such a condition may be regarded as being made up of short periods 
of quiescence followed by complete mixing, so that the final dissolved 
oxygen value at the end of any short period of quiescence becomes 
the initial, uniformly distributed value at the beginning of the next 
period. Since the rate of solution at the beginning of any period is 
strictly proportional to the deficit at that time, it follows that under 
the assumed conditions—and these are approximately the conditions 
met with—the rate of solution is proportional] to the existing satura- 
tion deficit. This is the law of the velocity of solution of solids in 
water, derived experimentally by Noyes and Whitney.” 


14 Noyes, A. A., and Whitney, W. R., Zeit, f. physik. Chem., 23, 1897, p. 689. 
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The physics of reaeration phenomena have been studied exten- 
sively by Adeney ** and his coworkers, who have developed the fol- 
lowing formula expressing the rate of reaeration of quiescent col- 
umns of water: 

w=(100—w,) (1- ett) 
where— 
w=amount of gas dissolved, expressed as percentage of saturation. 
W,= initial concentration. 
f=coefficient of escape of gas from the liquid per unit of area 
and volume. 
v=volume of liquid. 
a==area of surface. 
t=time of exposure. 

If this formula be differentiated with respect to time (t), it 

becomes— 


a = (100 — w,) -f(a/v)-e fe/v)t 


whence the rate of reaeration (SE) is shown to be directly pro- 


portional to the saturation deficit (100—w,). It is thus apparent 
that whether reaeration be viewed as a phenomenon of diffusion or 
one of “streaming,” both theories lead to the same conclusion as 
1egards the fundamental importance of the law of solution in deter- 
mining the rate of progression of the reaction. 

An experimental confirmation of this law as applied to the solu- 
tion of atmospheric oxygen by water was obtained by Dibdin™ from 
a series of about 150 tests made in connection with experiments on 
the condition of the water in the River Thames, in England. He 
exposed deaerated water in open vessels for periods of time rang- 
ing from 1 to 96 hours, determining the dissolved oxygen at in- 
termediate intervals. The results obtained were plotted in the 
form of a curve, which has been reproduced in a different, form 
in Figure No. 3, the ordinates being saturation deficit values plot- 
ted on a logarithmic scale and the abscissae, corresponding times 
from the starting point. Jt is noted that the plotted observations 
lie almost exactly along a straight line. Denoting the time as (t) 
and the saturation deficit as (D) we have, then, that— 


log D = —(ct + d) 


1% Adeney, W. E. Sci. Proc. Royal Dublin Soc., 1914. 

Adeney, W. E., and Becker, H. G. Philosophical Magazine, vol. 38, 1919, pp. 317-338. 
Philosophical Magazine, vol. 42, 1921, pp. 87-96. 

Adeney, W. E., Leonard, A. S. G., and Richardson, A. Sci. Proc. Dub. Roy. Soe., vol. 
17, 1922, pp. 19-28. 

16 Dibdin, W. J.. The Purification of Sewage and Water, 3d, ed., 1908, pp. 283-284 
(Diagram §8). 
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the constant (c) being the slope of the line and (d) being the in- 
tercept on the (D) axis. Differentiating this equation with re- 
spect to (t), we have— 
-S- cD | 
Since (Se) defines the rate of solution of oxygen in terms of 


saturation deficit, it is apparent that this rate is directly propor- 
tional to the saturation deficit (D), which is in accordance with 
the general law of solution above stated. 

From what has been previously stated regarding the relation of 
mechanical mixing to diffusion and “ streaming” as agencies in the 
reaeration of flowing bodies of water, it naturally follows that re- 
aeration of streams, while conditioned primarily by the rate of solu- 
tion of oxygen at the water surface, is modified by those factors 
which affect the rapidity and thoroughness with which the oxygen, 
once dissolved, is distributed throughout its depth. Thus it was 
found by Fair*’ that at a given saturation deficit value, the rate of 
reaeration of a quiescent body of water is greatly accelerated merely 
by the mixing resulting from induced wave action. In Sheffield, 
England, an activated sludge plant for sewage treatment has re- 
cently been developed in which the high rate of aeration required 
for oxidation of the sewage is obtained entirely by a mechanical 
mixing device. Examples of this kind, which are numerous, show 
conclusively the great influence exerted upon the rate of reaeration 
of flowing bodies of water by those forces of mixing and convection 
which may be summed up under the term “turbulence.” In run- 
ning streams the turbulence factor is highly variable and produces 
correspondingly varied effects upon reaeration rates. In a given 
stream stretch and under a given condition of flow, where the turbu- 
lence remains fairly constant, the rate of reaeration is a direct 
function of, and should be closely proportional to, the prevailing 
oxygen saturation deficit. Under conditions found in streams, there- 
fore, the operation of the law of solution is fundamental, and varia- 
tions in the rate of solution are governed largely by those physical 
characteristics of a given stream which cause different degrees of 
turbulence. 


THE OXYGEN BALANCE IN A STREAM 


The two opposing reactions, deoxygenation and reaeration, tend 
always to come to a condition of temporary equilibrium. If the 
water be nearly oxygen-saturated, and highly polluted, there is a 











17 Discussion of paper by R. H. Gould, ‘‘ The area of water surface as a controlling 
factor in the condition of polluted harbor waters,” Trans. Am. Soc. C. E., vol. 85, 1922, 
pp. 728-731. Cae 
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rapid rate of withdrawal of oxygen and a slow rate of replacement, 
resulting in a decrease in the available dissolved oxygen. As this 
value decreases, the rate of reoxygenation is correspondingly in- 
creased until it equals the rate of depletion, at which point the 
two reactions are for the moment in equilibrium and there is no 
change in the actual oxygen content. This equilibrium, however, 
Is ony momentary, for the decreasing oxygen demand of the or- 
ganic matter, resulting from its own oxidation, makes the rate 
of depletion correspondingly less and permits the gradual recovery 
of the dissolved oxygen up to its full saturation value. Under 
conditions of continuous or repeated pollution, however, an equilib- 
rium point may be reached at which the rate of reoxygenation is 
exactly equal to that of deoxygenation, and is so maintained. The 
importance of the reoxygenation factor itself and of its accurate 
experimental determination is therefore obvious. Upon this value, 
under any given conditions, depends the resultant oxygen condi- 
tion of the stream for a stated degree of pollution, or, conversely, 
the maximum amount of pollution compatible with any stated degree 
of oxygen depletion. | 

From the primary laws of oxidation of organic matter and of 
reoxygenation of a stream, the resultant general equation of stream 
condition may now be derived, expressed in terms of dissolved Oxy- 
gen. The application of the experimental data to this equation will 
then permit the derivation of the various constants which character- 
ize the stream in regard to its capacity to receive and dispose of 
sewage pollution. 

According to the argument which has been presented, the rate 
of change in the oxygen deficit is governed by two independent re- 
actions. First, the deficit increases at a rate which may be assumed #8 
to be proportional to the oxygen demand of the organic matter. 
Secondly, it decreases by reaeration, at a rate directly proportional 
to its own value. The two rates may be expressed in differential 
form thus: 


dL dD : 
Ae aE (see p. 5) - 
1D, 
and ae a 


in which— 
t = time of reaction, in days. 
Li= oxygen demand of the organic matter, expressed in terms of 
parts per million of oxygen. 
D = oxygen saturation deficit of the water, in parts per million. 


— 


18 Wor evidence supporting this assumption see later text, pp. 40 to 44. 
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i as rate of deoxygenation, in terms of oxygen saturation deficit. 
dD, CECE 

“dt ~ rate of reaeration, in terms of oxygen saturation deficit. 

K, = coefficient defining the rate of deoxygenation. 

K, = coefficient defining the rate of reaeration. 


The net rate of change in the oxygen deficit (D) at any time is 
equal to the difference between (or algebraic sum of) the two par- 







Fig No.4 
CHART SHOWING FUNDAMENTAL FACTORS INVOLVED IN DISSOLVE? OXYGEN CHAMGES 
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L* Oxygen demand of organic matter, 
R= Re-oxygenation in hme "er 
+ = Time of reaction. . BPdlL-dR*K,Ldt -K2Pdt 


K, = De-oxygenation coefficient 3o « KiL*K,0 
Ke Reaeration coefficient 


tial and opposing rates as defined above (see on No. 4) and may be 
pressed mathematically as follows: 


dD_ dD, , aD, 


deta he. 
whence— | = =K,L-—K,D 


which is a linear differential equation of the first order (sometimes 
called Leibnitz’s equation) having the general form— 

Ce te es 

fe 
The integrated equation derived from this differential equation 


defines the actual dissolved oxygen content of the water, expressed 
as the deficit below saturation, at any time, in terms of the other 
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factors. The full mathematical derivation of the integr ated equa- | 
tion is given in Appendix A, the resulting formula being— | 


D be cok (e—Kit — e- Kit) +D, e— Kat (1) 
in which— 
D,=initial dissolved oxygen saturation deficit of the water, in. | 
parts per million. 
D=saturation deficit, in parts per million, after time (t). 
L,= initial oxygen gins of the organic matter of the water, in 
parts per million. 
K,=coefficient defining the rate of deoxygenation. 
K,—coefficient defining the rate of reaeration. 
t=elapsed time, in days. 
e=base of Naperian or natural logarithms—2.71828. 
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Common logarithms may be used in obtaining values of (K,) 
and (K,), in which case the quantity (10) is substituted for (e). 
The quantities (D,), (D), (L,), and (4) may be expressed in terms 
of either parts per million or per cent of oxygen saturation; but 
should be stated invariably in the same terms. 

The type of curve defined by this formula is illustrated by curve 
A in Figure No. 5, which is based on an assumed simple case wherein 
the water is saturated with oxygen initially and all of the polluting 
matter enters a given stretch of the stream at or above its upper 
limit (Station A). In the curve shown, however, the ordinates are 
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in terms of the oxygen saturation value, which is equal to 100 per 
cent minus the deficit (D) as given by formula (1). The curve, 
which is typical of oxygen conditions frequently observed in streams 
below major points of pollution, has a distinct minimum point, 
where the rates of deoxygenation and of reaeration are momentarily 
equal. By differentiating equation (1) with respect to time, and 
placing the resulting expression equal to zero, the point of minimum 
oxygen content, in terms of time (t), is thus defined: 


eC xot = geal 1 Daa By | (2) 


It will be noted that in equation (2) it variable oxygen deficit, 
represented by (D) in equation (1), is absent. In order to obtain 
its value, equations (1) and (2) may be combined, giving the fol- 
lowing expression: 


Log D+K,t= Log| Ls = a (3) 
1 





This equation can be cleared of all terms in (t) and the maximum 
deficit expressed in terms only of the four constants, but the expres- 
sion obtained is unwieldy, and it is more convenient to solve first for 
(t) in equation (2) and then for (D) in equation (3). 

Similarly, the temperature function of this point may be obtained 
in a single equation, giving the time and maximum deficit as affected 
by variation in temperature, but this also is inconvenient for prac- 
tical use, and the result is obtained more directly by using (K,) and 
(K,) values independently determined for the temperature in ques- 
tion. The temperature effect upon the constant (K,) has already 
been given; that for the reaeration coefficient (K,), which remains to 
be determined experimentally, will be discussed later in connection 
with the application of the experimental data to the determination 
of that constant. 

The significance of the various terms in equation (1) is fairly 
obvious, excepting that of the reaeration coefficient (K,), the mean- 
ing of the deoxygenation constant (K,) having been previously 
discussed (p. 6). The reaeration coefficient (K,) is analogous to 
(K,) in that it defines a geometric rate of progression on a time 
basis; for example, if the value of (K,) be such that 20 per cent of 
the existing saturation deficit is satisfied by reaeration in the first 
unit of time, then 20 per cent of the remaining deficit will be satis- 
fied in the second unit of time, and so on. It differs from (K,), 
however, in not being a constant for a given temperature, as is true 
(or approximately so) of (K,). It has already been noted that the 
rate of reaeration of a body of water is modified to a large extent 
by its degree of turbulence, other things being equal. In flowing 
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streams, where turbulence undergoes wide variations according to 
velocity of flow, character, and slope of the channel, and other physi- 
cal factors, correspondingly marked differences are to be expected 
in rates of reaeration and hence in values of (K,), observed in differ- 
ent stretches of the same river, or even in the same stream stretch 
under varying flow conditions. On the other hand, the reaeration 
rate should be closely related to and governed by those conditions 
which influence turbulence of flow; hence values of (K,), as actu- 
ally determined in a given river stretch, should bear a close relation 
to the several measurable factors of physical condition which cause 
varying degrees of turbulence. With a given type of channel or 
flow condition, values of the coefficient (K,) should be well defined 
and characteristic for that type or condition. 

The practical significance of the reaeration coefficient in problems 
of stream pollution is therefore twofold. First, there is the local ap- 
plication, in which a set of (K,) values, once determined for a par- 
ticular stream, may be utilized to calculate its capacity for reaera- 
tion under any assumed conditions of future pollution, using for 
this purpose the three formulas that have been developed above. (See 
pp. 18-19.) Secondly, there is the more general application, wherein 
correlations of values of (K,) with certain measurable factors of 
physical stream condition may be employed to estimate the reaera- 
tion capacities of other streams for which these factors are known but 
in which no direct. measurements have been made of the coefficient. 
(See p. 64.) Obviously the point of departure for a study of this 
kind is the direct measurement of the value of (K,) in the stream. 
The methods which have been employed in the derivation of reaera- 
tion coefficients for a number of stretches of the Ohio River will 
therefore be described in some detail. 

Referring now to equation (1), if the problem were that of solv- 
ing for an unknown oxygen content (D), with all of the other fac- 
tors known, the matter would be one of simple substitution. ‘This 
would be the procedure in practice, the constants being given and it 
being required to determine the residual dissolved oxygen after any 
time and under given or assumed conditions of pollution. In experi- 
mentally measuring the reaeration coefficient (K,), however, all of 
the other terms, including (D), must first be determined. For con- 
venience in reference, these terms will be listed again, as follows: 

K,=the deoxygenation constant. 

L,=the initial total oxygen demand of the stream water in parts 
per million of oxygen. 

D,=—the initial dissolved oxygen content of the stream water in 
parts per million, expressed as saturation deficit. 

D=the dissolved oxygen content of the stream, in parts per mil- 
lion, after a time “t” from the initial point. 
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t—the time of flow, in days, from an oxygen content of (D,) to 
one of (D). 

The determination of (K,) and the temperature relations involved 
already have been discussed. In the present case values have been 
assigned to (K,) by correcting its value, 0.100, at 20° C. to its 
various equivalents at different observed stream temperatures, 
using for this purpose the curve shown in Figure No. 2 or the cor- 
rection formula given on page 7, with the thermal coefficient (9) 
taken as equal to 1.047. The determinations of (D,) and (D), and 
of the time (t) in a given stream stretch, require the selection of two 
sampling stations, A and B, one located at the upstream and the 
other at the downstream end of the stretch, and direct observation 
of the dissolved oxygen content and the water temperature at these 
two stations, together with measurements or computations of the 
mean time of flow between them at the various river stages at which 
the dissolved oxygen content is observed. 

The assignment of a proper value to the initial oxygen demand 
(L.), however, presents an extremely difficult problem, owing largely 
to the uncertainties involved in determining, under conditions found 
in natural streams, a quantity which may be taken as being a rep- 
resentative one for a particular river stretch under observation. It 
is proposed, therefore, to discuss this question in considerable detail. 

If the river stretch in question were entirely free from inflowing 
pollution or dilution at points intermediate between two given 
sampling stations, or if the amounts and points of entry of such 
inflow were definitely known, the various factors concerned in the 
problem would be determinate and its solution would be compara- 
tively simple. In the first case, the only change in the initially 
observed oxygen demand, occurring between the upper Station A 
and the lower Station B, would be a progressive and orderly decrease 
due to oxidation, following a course which, for all practical pur- 
poses, could safely be assumed to be similar to the logarithmic time 
function curve defined by line A in Figure No. 6a. The initially 
observed oxygen demand, denoted as (La) in the chart, therefore 
would become the value of (l.) in formula (1). In the second 
case, where known increments of pollution or dilution entered the 
river stretch at recognized points between Stations A and B, the 
position of the line of residual oxygen demand would be altered at 
each one of these points, as shown in Figure No. 6a, but its slope, 
defining the rate of oxidation, would remain constant for the partic- 
ular river temperature condition prevailing. In this latter instance, 
a new value of (Lz) would be computed at each point of entry of 
known inflow, the new ordinate of the curve lying above or below 
line A according to whether the effect of the added increment of 
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overflow at that point were one of added pucenare or added dilution, 
respectively. 

The condition actually encountered in a large majority of cases, 
however, is far from being the simple one above described, owing to 
the fact that the density of oxidizable matter carried by natural 
streams is subject to frequent disturbances due to the effects of in- 
flow, sedimentation, and channel scouring. The combined influence 
of these disturbing factors may bring about marked changes in the 
oxygen demand of a stream wholly unrelated to those which may be 
caused by the orderly processes of oxidation; thereby masking the 
latter, in so far as they may be observable, vid rendering uncertain 


Fic.No.6a 
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Line A+ Oxidation undisturbed by inflow. 
K,: Rate constant at stream temperature. 
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any attempted determination of a representative oxygen demand 
figure for a particular river stretch from observations at fixed sam- 
pling stations. In the absence of specific information as to the 
amounts and points of application of the various disturbing factors 
noted, the only clue as to their combined influence is afforded by the 
extent to which the oxygen demand (Li), observed at station B 
(fig. No. 6a), deviates from the residual component (L’,) of the 
initially observed oxygen demand (Lia) remaining at station B, 
assuming (La) to have been satisfied along line A in accordance 
with the law of oxidation stated on page 5, and making due allow- 
ance for any deviation due to errors of observation or extrapola- 
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tion.1® The extent of divergence between (Lg) and (L’,) gives no 
indication, however, as to the total amount of intermediate disturb- 
ance or the manner in which it may be distributed between A and B. 

Under these circumstances, any value assigned to (L,) for substi- 
tution in formula (1) must be an hypothetical one, based on some 
method of combining the results of observations at stations A and B 
in such a way as to give, as nearly as practicable, a measure of the 
total amount of oxygen demand satisfied between the two stations. 
With the foregoing limitations, it is fairly evident, however, that no 
single value of (L,), assigned according to a standardized procedure, 
can exactly satisfy this requirement. Any procedure adopted, there- 
fore, must be a compromise at best, designed to secure a fair ap- 
proximation to the result sought. 

There are three possible methods by which a value of (L,) may 
be derived from observations of oxygen demand at A and B. Thus, 
{L.) may be based on: 

1. A mean * of the two observations (ia) and (Leg). 

2. The observed value either of (ia) or of (Lz). 

3. A mean? of the observed (Lg) and the residual component 
(L’,) of the observed (Lia) remaining at station B. 

Method (1) obviously is illogical in that it assumes a straight-line 
change in oxygen demand from a value of (La) at station A to one 
of (Lg) at station B and fails to take any account of the expected 
progressive decrease in (lia) between the two points, due to oxida- 
tion. The sources of error which method (2) entails may best be 
appreciated by referring to Figure No. 6b, in which line A and line 
B illustrate the assumptions involved in adopting, respectively, an 
observed value of (La) or one of (Lg) as the hypothetical value of 
(L.). If (La) is used for the purpose, any divergence between (Ls) 
and (L’,) at station B is considered negligible, and a condition of 
zero disturbance is thus assumed. As previously noted, such a con- 
dition rarely, if ever, exists in natural streams. If (Ls) is em- 
ployed, the amount of divergence between (Lg) and (L’,) is as- 
sumed to be due to an initial oxygen demand equivalent to the quan- 








19 Errors of observation, or those involved in extrapolating the observed quantity (La) to its residual 
value (L’a), at Station B, may cause a divergence of (Ls) from (L’a) entirely dissociated from any devi- 
ution due to intermediate disturbance. The effect of observational errors should, in general, be so small 
is to be inappreciable. The effect of extrapolation errors, by which is meant principally those due to 
rrors in the constants employed in calculating the value of (L’a), might be sufficiently great, where the 
ime factor were long, to cause an appreciable difference between this quantity, as calculated, and the 
ybserved value (Lp), owing to the fact that such errors become enlarged, in comparison with (Lx), when 
xtrapolated over long periods of time. Any slight error in the assumed value of (Kj), for example, or in 
ts temperature correction factor (6), might produce an effect of this kind if magnified over a long time 
yeriod. For this reason, observations in river stretches of considerable length, where the period of flow 
xceeds several days, should be avoided where possible and their results interpreted with due regard for 
his source of error. 

2 The term ‘“ mean” is used here as signifying a simple arithmetical average. A 
weighted mean would hardly be justified in either case, owing to the absence of any 
lefinite basis for establishing relative weights. 
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tity (L’s) in the chart, representing the maximum amount of dis- 
turbance which could possibly account for the divergence observed 
in a given case. The use of either (La) or (Lg) as a basis for (La) 
represents, therefore, the two most extreme assumptions possible re- 
garding the amount of intermediate disturbance. In the one case, 
the derived value of (L.) evidently would be too low and in the 
other too high. Sale 

Method (8) represents a compromise between the two alternative 
procedures described under method (2), the primary assumption here 
involved being that the total. amount of disturbance between Sta- 
tions A and B is equal to an amount which, if concentrated at a 
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} CHART SHOWING METHOD OF OBTAINING HYPOTHETICAL 
VALUE OF Oo) BY AVERAGING OBSERVATIONS TAKEN AT 
AYANDP (B)AND BROUGHT TO A COMMON BASIS 

FOR COMPARISON AT STATION (6) 





STATIONS 


Note 

Line A-Oxygen demand curve based on observation at Sta.A 

Line B: Oxygen demand curve based on observation at Sta.B 

Lime C*Oxygen demand curve based on mean quantity (Lm) 
Slope of each line determined by assumed oxidation 

rate constant (K,), at stream temperature. 











SiocHEMicAL OXYGEN [DrEMmManp, P.P.M. 





Time OF FLOW From STA.A TO STA.B 


point immediately below Station A, would account for exactly one- 
half the observed divergence between (Ls) and (L’,4) at Station B. 
The effect of this assumption is illustrated in Figure No. 6b by line 
C, the value (L,) being an arithmetical mean of (Lg) and (L’,), or: 
Tet de a 
2 

It will be noted that the total amount of oxygen demand satisfied 

between Stations A and B along line C is equal to the arithmetical 


mean of the total amounts satisfied, respectively, along lines A and 
B; that is: 


La= 





| —L’ L’,—L, 
Clas — Lim) aes (L, a) 4 ( p) 
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Method (3) thus involves a further assumption that the total 
amount of oxygen demand satisfied between Stations A and B is 
equal to the mean of the amounts which would be satisfied under the 
maximum and minimum possible conditions of disturbance, respec- 
tively. In the absence of any specific data as to the actual amount 
or distribution of such disturbance, an assumption of this nature 
would appear to be a far more reasonable one than those which are 
involved in methods (1) and (2). The calculation of (L,) by this 
method is simple, its value being given, in terms of the observed 
quantities (La) and (Ls), by the relation: 


Cig lhe ta. ay ee a 
Se 


The methods above described have reference solely to the deriva- 
tion of a value of (L,) for substitution into formula (1) which, it 
will be recalled, is based on a primary assumption of an initial oxy- 
gen demand, subsequently undisturbed in its orderly rate of decrease 
through oxidation. It may be argued that procedures based on such 
an assumption can hardly be apphed, logically, to conditions of nat- 
ural streams in that they fail to take direct account of the manner in 
which the disturbing influences causing a deviation of (Lg) from 
(L’,) are distributed along the stream between two given observation 
stations. A contention of this kind unquestionably is sound theo- 
retically, but when an attempt is made to develop a procedure based 
on some definite primary assumption as to the distribution of dis- 
turbing infiuences along a stream, two difficulties present them- 
selves. First is a question as to the kind of distribution which shall 
be assumed in a given case, with little or no information available 
on which to base a particular assumption. Second is the obvious 
fact that for each different assumption made a new and character- 
istic equation of the type of formula (1) must be developed, starting 
with a new differential equation of condition, and undergoing the 
various steps of integration similar to those outlined in Appendix 
A with respect to formula (1). 

The mathematical complexities which such a procedure involves 
are well illustrated by the example, given in Appendix B, of the 
derivation of a new resultant oxygen formula based on the most 
elementary assumption possible regarding distributed disturbance, 
namely, a uniform distribution of inflow along the stream between 
Stations A and B. The equation thus derived is as follows: 


Kip Ror St)! (Ties eK iy, aye: ea: e~Ket)] 
K, (K,—K,) (—e7**) 

the symbols employed being the same as in formula (1), except that 

(La) and (Ig) represent the observed oxygen demand at Stations 

A and B, respectively. 


De 


+ D,e~K# 
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Efforts to develop other formulae of the same type, based on other 
assumptions as to the distribution of inflow, have led to integrated 
equations even more complex than the one above stated. For pur- 
poses of practical application such formulae are so complicated as 
to be virtually unworkable in most cases. In the absence of specific 
information as to the manner in which disturbing influence may be 
distributed in a given instance, the use of a particular formula, 
based on a single arbitrarily assumed distribution, would be hazard- 
ous, moreover, in view of the extreme variability and high degree of 
uncertainty with which such distributions occur in natural streams. 
With the admitted theoretical advantages of these procedures over 
the more simple one embodied in formula (1) their practical disad- 
vantages are of such a nature as to render their use inadvisable for 
purposes of working calculations. 

For use in connection with the analyses of the Ohio River data, 
tobe discussed in the text which follows, the procedure finally 
adopted as being the most suitable for the purpose at hand was the 
employment of formula (1), deriving values of (L.) by averaging. 
those of (up) and (Lye-K,t), as described under method (3) on 
pages 23-24. By this comparatively simple procedure values of (Li) 
were calculated from the base data given in Table No. 1 for each 
stretch of the river for which observations were available. 

An inspection of formula (1) indicates that a direct solution of it 
for the value of the reaeration coefficient (K,) involves mathematical 
difficulties, since this term appears both as a coefficient and as an 
exponent. While it is possible that a convenient solution of this 
equation by means of a nomographic chart might be devised, a farrly* 
simple method in practice has consisted of the indirect procedure 
of assuming values of (K,), solving for the corresponding (D) 
values; plotting these on cross-section paper with (K,) values as 
ordinates and (D) values as abscissae, and from a smooth curve 
drawn through the points, selecting the value of (K,) correspond- 
ing to the known (D) value. The substitution of the value of (K,), 
thus obtained, into formula (1) will insure that the proper figure 
has been chosen from the plot, thus checking the work. This method 
of procedure is comparatively simple in routine work if the so-called 
“log-log” slide rule is used, permitting the convenient solution of 
power functions. It is rarely necessary to obtain more than three 
points on the smooth curve to obtain the interpolated value of (K,) 
sought. 

Employing the foregoing method, values of the reaeration co- 
efficient may be readily computed for any river stretch or any con- 
dition of flow or season for which the necessary laboratory and 
hydrometric data are available. From calculations of this kind a 
set of (K,) values will be obtained for each particular river stretch 
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studied, the values varying according to the different conditions of 
temperature and flow encountered during the period of observations. 

The coefficients thus obtained, representing in each case the co- 
incident stream conditions, are modified by three major factors, 
namely, temperature, stream depth, and turbulence. The tempera- 
ture influence is, strictly speaking, a compound one, depending upon 
the resultant effect of temperature variations upon the rate of solu- 
tion of oxygen at the water surface and upon its velocity of dif- 
fusion in the lower strata. These two influences are opposed to 
each other in the sense that the rate of solution diminishes with ris- 
ing temperature, while the rate of diffusion increases. 

The effect of temperature variations upon the rate of oxygen 
solution at the surface is governed by the relative amounts of satura- 
tion deficit represented by a particular dissolved oxygen content 
when referred to saturation values at different temperatures. If, 
for example, the oxygen content be 5 parts per million, its satura- 
tion deficit at 20° C. is 4.2 parts and at 30° C. is 2.6 parts, the rela- 
tive rates of solution at the two temperatures being defined by the 
proportion 4.2:2.6. Where the oxygen content is expressed, how- 
ever, in terms of saturation deficit, as is the case in formulas (1) 
and (3), this temperature relation is taken account of automatically, 
and the diffusion relationship is the governing influence. 

In their New York Harbor studies, Black and Phelps?! derived 
experimentally a curve showing the effects of temperature varia- 
tions upon rates of diffusion of atmospheric oxygen in water. In 
Figure No. 7 this curve, slightly modified for temperatures below 
10° C., is reproduced in such a way as to show relative rates 
of diffusion at various temperatures with reference to the rate 
of 20° C. From this curve it appears that the rate of diffusion 
at 20° C. is doubled at a temperature of about 27.2° C., trebled 
at 31.4° C., and halved at 10.6° C. It is evident that the tem- 
perature of the water, even within ranges ordinarily observed in 
streams, exerts a great influence upon the rate of diffusion of oxygen 
in water and hence upon the rate of reaeration of streams. By 
means of a curve similar to that of Figure No. 7, values of the 
reaeration coefficients as determined at various prevailing stream 
temperatures may be reduced to equivalent values at a standard 
temperature of 20° C. before comparing them with reference to 
other stream conditions. 

The influence of stream depth upon the rate of reaeration is con- 
trolled in part by the relation of depth to volume of flow (which, 
in turn, governs the concentration, in the stream, of oxygen derived 


21 Black, Col. W. M., and Phelphs, E. B., The Discharge of Sewage into New York 
Harbor. Report made to the Board of Estimate and Apportionment, New York City, 
1911. 
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from reaeration) and probably also in part by the relation of depth 
to the rate of diffusion. It happens that both the rate of diffusion 
and volume of flow are power functions of the depth, approximat-’ 
ing its square. A glance at the equivalent of (K) in the diffusion 





SAH et ‘abr : 
formula on page 11, which is given by the quantity (4m) will 


show, for example, that the rate of diffusion is inversely propor- 
tional to the square of the depth (L). As regards the volume of 
flow, it is equal to the velocity times the area of flow, the relation 
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CHART SHOWING RELATION OF VALUES OF REAERATION COEFFICIENT 
K, AT VARIOUS TEMPERATURES TO ITS VALUE AT 20°C. 


Derived from Black & Phelps, Report on Discharge of Sewage 
into New York Harbor, Page*S6. Curve slightly ‘modified 
for temperatures below 10°C 


RATIO OF VALUE OF Kz To mS VALUE aT ZO°C 





TEMPERATURE OF WATER IN °C, 


being expressed by the simple hydraulic formula: (Q=AV), in 
which (Q) denotes the volume of flow, (A) the area, and (V) the 
velocity. In most streams having a fairly large ratio of width to 
depth, the area of flow (A) and the velocity (V) are each of them. 
very nearly proportional to the depth (H), so that the relation to 
volume of flow may be written as being roughly: (Q=cH?). The 
amount of reaeration per unit of time, or its rate, when measured in 
terms of oxygen concentration, should be inversely proportional to 
the quantity of water throughout which the oxygen is distributed; 
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that is, to the stream discharge (Q). From the relation just shown, 
this rate should therefore be roughly proportional directly to the 
square of the depth (HH). 

As to the relation between reaeration and turbulence, no well- 
defined theoretical basis exists for judging its nature, the reason be- 
ing that the turbulence of a stream is an abstract function dependent 
upon and capable of expression only in terms of some other more 
concrete stream characteristic. For a given stream type, some 
definite relationship might reasonably be expected to exist between 
the turbulence of a stream and its velocity of flow, but it is evident 
that different types of streams should reveal quite different velocity- 
turbulence relations. For example, in a stream possessing an even, 
smooth, deep channel, with flat slopes and either straight or gradu- 

ally changing direction of flow, the effect of velocity upon turbu- 
lence is conceivably much less than in a watercourse having a shal- 
low, rough channel, with steep’ slopes and sharp changes in direc- 
tion. Between the two extremes lie numerous well-defined stream 
types; and a large river such as the Ohio is likely to exhibit through- 
out its course several of these types. 

No very definite theoretical relation exists between turbulence and 
velocity, even assuming other physical conditions constant. In a 
general way, turbulence is the result of frictional resistance to flow, 
and under uniform physical conditions might be expected to be a 
power function of the velocity of the form: 
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the constants c and n defining the stream type as regards the fixed 
physical conditions, such as slope, character of bottom, depth, shape, 
and direction of channel, etc. This subject can therefore be dealt 
with only empirically. 


EXPERIMENTAL 


The quantitative working theory of the process of oxygen stream 
purification which has been developed and expressed in the form of 
certain algebraic equations is capable of two distinct uses. With 
all the constants known, it makes possible the computation of ac- 
tual stream conditions and may even be employed for the determi- 
nation of future conditions, with an assumed increase in con- 
_ tributing population or an assumed degree of sewage purification. 
Information is lacking, however, concerning the reaeration or re- 
oxygenation coefficient (KX,) in the formulas, excepting from a 
purely theoretical standpoint. 

It is possible, however, to utilize these formulas for the experi- 
mental determinations of this term, by actually measuring the 

999—25——_3 
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more readily determinable present stream conditions. To this end 
analytical data have been obtained in the form of dissolved oxygen 
and oxygen demand values, and these data will be employed for the 
determination of the reaeration coefficient. It will then be possible 
to study this coefficient in connection with the physical stream con- 
ditions which affect it and to search for any general relationships 
which would make the determinations of reaeration: a matter of 
computation for known physical conditions. | 

It should be noted that the other data necessary for the solution 
cf the condition equation are either known or capable of ready deter- 
mination in any particular case. 


PRESENTATION OF BASE DATA | 


The scope of the observations constituting the basis of the study 
embraced daily or thrice-weekly determinations of dissolved oxygen 
and biochemical oxygen demand of the Ohio River at selected stations 
throughout its entire length, from Pittsburgh to below Paducah, 
over a period extending from May 1 to October 15, 1914, supple- 
mented by similar observations in a stretch of the river extending 
from above Cincinnati to below Louisville, over a period from Octo- 
ber 16, 1914, to April 30, 1915. The data for the latter named stretch 
thus embraced a full year’s cycle of stream conditions. | 

The base data derived from these observations were first reduced 
to terms of monthly average figures and summarized as shown in 
Table No. 1, which. also includes certain supplementary data 
employed in the analysis of the laboratory results, such as mean river 
water temperatures and times of flow between the various sampling 
stations (columns 2, 3, and 4). The sampling station notations as 
given in column 1 require explanation. The location of all stations 
was referred to a point at the junction of the Allegheny and Monon- 
gahela Rivers at Pittsburgh, and each station on the Ohio River was 
given a number corresponding to its location in miles below the ref- 
erence point. (See Map A, showing the location of the various sta- 
tions; also Public Health Bulletin No. 1438, pp. 98-104, for a descrip- 
tion of the stations.) Thus the station designated as “ Ohio No, 23” 
was located 23 miles downstream from the reference point. The tribu- 
tary stations, designated by name, were located in all except two 
cases practically at the tributary mouth. The two exceptions were 
the Allegheny and Monongahela stations, which were located at 
distances of 7 and 12 miles, respectively, above their junction at 
Pittsburgh. 

In the analysis of the data which follows, 11 stretches of the river 
have been selected for reaeration study, the various factors involved 
in the calculations having been transcribed from Table No. 1 and 
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arranged as shown in Table No. 2. The last two columns in these 
tables contain values of the reaeration coefficient (K,) as determined 
for the actual stream condition by the method outlined on pages 
26-27 and as converted to a standard temperature basis of 20° C., 
using the curve of Figure 7 for the purpose. 
is given a summary of the values of (K,) thus standardized and 
brought together from Table No. 2 for reference and comparison. 


In Table No. 3 


TasLeE No. 1.—Average results, by months, of laboratory determinations of 
dissolved oxygen and biochemical orygen demand at various sampling sta- 
tions i the Ohio River and at the mouths of designated tributaries, together 
with collateral data concerning mean river water temperatures and times 
of flow between successive stations 
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TaBLE No, 1.—Average results, by months, ete-—Continued 
MONTH OF JUNE, 1914—Continued 





Mean time of 






























































DODO PNWDOONWRDOrHO 





























ae F , Total oxygen 
foe: days: Initial dissolved oxygen rae darnkiod 
during 
a $4). |b ee 
i Soa . hours’ | Quan- 
‘ incu- tity 
; ; tem- é ia 
Sampling station pera- Ban a oe bation hath 
a- er + ration a arts 
og Pitts- | tion ened cent | deficit | 20° C. Parts sie mil- 
- | burgh | next. |P lion | Satu- | (parts | (parts | hic lion X 
above ration per mil-/per mil-- 7°" | ¢hou- 
lion) lion) | ; sand 
late second- 
; | feet) 
(1) (2) (3) (4) (5) i) a ls (8) *} > 8) (10) 
@InIGRINIGS BAG 2 ees a oe Dba 17. 82-1 LOL09 7.42 | 89.4 88 542 2. 04 36. 
SClOtOse noe Sete ees 24.6.) 17. 52 . 20 Tad). OOeG: | 69 1. 03 5. 00 5. 
Onto, No, 401 SS et 26.2 | 20.24) © 2727s 8B iy 96.3 21 60-53: 48 | ae 
TICKING se aon wee 27.6 | 20. 54 . 30 6.88 | 86.2 1 SEO . 98 4. 76 4, 
ORIGIN. 402 ar. tates. oe 26..2°| 20. 91 .387 6. 35 77.5 1. 84 . 90 4.37 91, 
Ohio;.No; 482. see 2... 3 25. 8 21. 39 . 48 6.7 82.3 1. 46 .91 4, 42 92. 
OHO, NO. 488-5. Fes SSIS a6 6 | a ae o20 6. 45 77.2 1. 838 . 74 3. 60 75. 
Winwrawie: ey 8 See 25.6 | 21.69 . 05 7. 90 95. 4 38 1 oO) 7. 28 9. 
One: NO. 4920 2 BIR 25. 0 PA br Ae . 08 6. 48 77.3 1. 90 . 90 4. 37 98. 
Obit No, 598.22: BOeee wae bass, sereoo 6. 06 7. 51 90. 2 82 . 95 4.61 132. 
Cag NG, Bite soe. Las 25.9 | 28.78 . 95 7.3 89. 6 86 . 82 3. 98 114. 
ORIGIN 0. 904. 22.2. aL ae 26.91? 412 62. || 138209 7. 24 89. 6 84 . 40 1. 94 on, 
@uimpber land... 22 eu ee 29.4 | 41.96 . 34 5. 56 | 88. 3 2.15 . 00 . 00 : 
OIG? NOs 920. 2.4 Sa et 26.9 | 42.34 . 88 7.44 | 92.1 64 eon L70 86. 
"MemiessOG= =. 2. wees oe 29.6 | 42.40 . 06 6.56 | 85.3 1.13 . 20 1. 26 PAS 
ORIOFNo,9838_.2. = 2 Nits alk eae 28.0 | 42.95 . 65 1220 90. 9 oe ie 3. 74 255. 
$e i 
MONTH OF JULY, 1914 
We Sh. Ba Se ua ee es Se ee a ee ee ee ee 
t | | 
Allegheny; NoOw7 sib. 22. ss Dde Gels ee pee a oe ees 7. 36 8% Ba 71.02 0. 33 1. 60 65 
Monongahela, No. 12_._.---- aie | a ee a 6. 58 78. 6 1. 78 43 2. 09 4.0 
Ono; N@r8e OF ae ost ae Piet 0. 68 0. 68 4. 36 52, 1 4. 00 44 2.14 DH5 ; 
Obie, Noell. 4 sa... <5 eee 24. 0 2. 49 1.81 7. 21 84. 5 1.32 .41 1.99 10.7 
Oi, Now] 922. seks. Sees 24, 2 3. 85 1. 36 7. 46 87.8 1. 04 ~ 44 2.14 de5 
Onio; Now23....2.0e. =) fae 24. 5 4, 45 . 60 7, 21 85. 3 1. 24 . 48 2,33 12.6 
HOA VG? iy Gee. eae ee eee 23. 8 4,74 29 742° 86. 7 1.14 . 64 3. 11 1.4 
Ginko: IN GsbGe 3. See Dee 24. 6 7.45 2.71 7.50 |. 88.9 . 94 . 24 1.16 RT 
ONTO NORT 4.3 eas 5 ee 25.0 | ‘9..09 1.73 ‘7,84 93. 5 54 . 39 1.90 12.6 
(rho; INOMSSe2 ss Beef. ak ee 25.0 | 10.20 1.01 7. 86 93.8 | 52 . 49 2. 38 15.8 
Ontos NOsO 72. 40852255524 24.8 11. 20 1. 00 7.74 92. 1 67 .33 1. 60 11.6 
Onio, NowlO4....8- 2.) Bt. 24 Sah 12; dl 91 7.49 89. 0 . 92 AS! P3351: 10.0 
Odio, Now840..- 230i4.2.4. 424 - 26:74. 24°90 | 12.79 6, 83 84,21 - 2,28 . 38 1,85 30.3 
SRC te a sa Se 2h Gul. 25.40 20 7.57 89. 6 | 71 . 90 4, 37 3.2 
Olio, Noid@l 2a. 2.5 the 27.1 27. 86 2. 76 7.70 95. 6 35 . 40 1.94 36.3 
AGRO oe Go). Bn 2 dhe 29.3 | 28.21 35 6. 63 85. 7 1.10 . 83 4. 03 .4 
OD10; NOAA. 222. Gin. ccke Be 27.3 | 28. 64 43 6. 00 74.9 2. 02 . 81 3. 94 74.8 
Ohio, NO#482:. - 64-5. .22 ee 26. 8 29. 20 56 6. 76 83. 5 1. 34 . 81 3. 93 74.6 
ODO; AN OSAS8 S22 oe eee ee 27.2 | 29.46 26 6. 62 82. 3 1. 42 . 68 3. 30 62. 7 
CURTO, NG 4G2: owe So 5 26.5 | 29. 60 14 6. 43 77.8 yal . 66 3, 20 64. 0 
Ohio; No, 1508.32 25.- 32S 27.0 | 36.86 7. 26 7.70 95. 5 GA 1.01 4, 90 113.0 
ONO INONGH 3 222k. oak eS 27.2 | 38.06 1. 20 7. 45 92. 6 . 59 . 84 4. 08 94, 4 
Ohio IN G#6ROc: «SS 2- se se 27.2 | 38.38 32 7. 06 87.8 . 98 . 98 4. 76 109. 8 
Ohio, NGow904.. .ce<. 252 sts 28.3 | 53.19} 14.81 7.61 95.5 . 36 382 1. 55 47.7 
GCumberland:...s222-. 22 lies 28.9 | 53.64 45 6. 28 80. 8 1. 50 . 26 1. 26 10.1 
OIG MNO. O20 bs 225. 3. sons 28.9 | 54.05 41 7. 34 94.4 . 44 .21 1. 04 40.3 
MGW OSSOO/ Ao a Socket cuceas ee 29.7 | 54.11 06 6. 47 84. 3 1. 20 . 18 Wel! 16.7 
OHio eNOn088 2235-655 secant 29.0 | 54.69 . 58 6. 90 88. 8 87 . 28 1. 36 78.8 
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TABLE No. 1.—Average results, by months, ete —Continued 








MONTH OF AUQUST, 1914 
























































Mean time of 
: Total oxygen 
sa — Initial dissolved oxygen meee demand 
during 
Mean ie home? Quan- 
water ineu- tity 
5 " tem- A i 
Sampling station cia Satu- | bation units 
Sa Sta- Rats Per | ration at Parts (parts 
og” | Pitts- | tion jo. mil-| cent | deficit | 20°C. per mil- 
‘ | burgh | next Oe on satu- | (parts | (parts | Pijon| lion X 
above ration |per mil-|per mil- thou- 
lion) | lion) sand 
seconde 
feet) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 
AllechenveNe: 7223-84-82 25; Osan. 2 ees Se HES 82. 4 1. 52 0.39 1.89 3.1 
Monongahela No. 12......--- P21 a MS eae ty Eee 6. 40 76.0 2.01 en 3. 45 7.1 
GistouNG: ome: | Ve 5 he 25.0 1.00 1.00 2. 85 34. 0 5. 53 .61 2. 96 11.0 
@ito-Nob best. ae. Pk 24. 5 3.73 2:73: 6. 50 76.9 1. 95 OL 2. 48 9g. 2 
Bene: Nee 19s. 5 Bo... We 24.4 je7t 1.98 6. 87 81.1 1. 60 woe 2. 52 9.4 
Whio NO: 2o-)..235. 33... 5 BS 24.0 6. 59 . 88 6. 65 Vie) 1. 88 . 38 1. 85 6. 9 
IBGavete wee Re cet 23.0 7.01 . 42 7. 40 85. 2 1. 28 702 2.52 9 
ROUND: 6522... 45... 2 we 24. 2 11. 49 4. 48 7.75 91.2 Bhs) . 28 1. 36 6.3 
Obio NOs Wisden Se 3 24,6 13. 81 2.47 7.95 94,2 .49 .12 . 58 257 
@imo Not S82... Feo. 24. 4 15. 47 1. 51 7.76 91.7 71 . 50 2. 43 11.3 
Olno Noo Ora. 2 oS 24, 2 16. 92 1.45 7.90 93. 0 . 50 7Oo 2. 67 12.4 
OmMo Noe 1044 =. 3 Sh. A 23.9 18. 26 1. 34 7. 59 88.9 . 95 . 35 1.70 7.9 
Ohio Nor 340! 2 4-38. 2 20.94 34. 88 16. 62 7.09 85. 8 1.18 =o 1.80 21. 2 
SelotvOsa: Be eo 5 Be gs 2 A 24.8 30; Ee . 24 (et 87.3: 1.04 .83 4.03 3.9 
Ohio No. A6ls 432k 26.4 | 38.34 3. 22 (evs: 89.1 | . 89 .49 |. 2.38 33. 0 
Mielainigey ee. Gales 22s US 27.2 | 38.64 . 30 6. 58 81.8! 1.46 1.01 4,90 3.6 
Olio Not 4752.23.28 30 1 SR 26.3 39. 32 . 68 4.97 60.8 3. 20 79 3. 84 58. 4 
Ohio No: 482> .. s- eho: BA 25.8 40. 01 . 69 6. 39 Tb 15-86 wt4 3. 59 53. 0 
Obio-No. 4885... 5 Sto. ie 25.9 40. 30 . 29 6. 24 75.8 | 2.00 | . 65 3.15 47.9 
Olio: No. 4922.5 2 Bee Ee 25. 0 40. 45 . 15 6. 23 (EROS IE eats . 66 3. 20 52.8 
Ohio No; 698s = 2-165. oe E 26.0 48. 98 8. 67 7.81 96. 6 . 41 “30 1. 60 29.3 
Oino No: 6lis 2. 22. ae - 27.0 | 50.42 1. 44 7.79 95. 0 28 . 38 1.85 33.8 
Ohio No. 6492... 5-222... 2 St Zhai 50. 80 . 38 Go83 91.1 72 - 60 2.91 53. 2 
Obie Noy 0045 - 25 Ee... ER Dio 67. 54 16. 74 7.18 89.9 . 81 . 16 . 78 20. 6 
@umberlands..--i2s2...24:2 28.0 67, 88 . b4 6. 28 79.3 1. 64 | 18 . 87 6.0 
@hio: No: 920: =. Pf... ae 27.9 68. 37 49 7.78 98. 2 ai Bye . 40 1.94 64. 5 
‘PEMNCSSOO eo. 5 265... 2h 28.4 | 68. 43 . 06 6. 95 88. 4 91 24 alisra te 19.5 
Ointo No. 983: . 2.20.2. b SL 27.9 69. 00 Vow 7.09 89. 4 . 84 | . 30 1. 46 73. 0 
MONTH OF SEPTEMBER, 1914 
Alteghenyi NO. 7.-25.2....28. 19.0 | Rapa as 5a a Ga 8. 25 88. 2 1.10 0. 46 2. 24 5.0 
Monongahela No. 12..._.___- QAO Rie: = ee 6. 54 71.3 2. 45 ~ 55 2a6u 3.0 
OlaOmING. Bess. . 25 28825 k Te 2150 1,05 1. 05 Pie 80. 2 6. 28 . 59 2. 86 9.6 
Wivtor NG. 625 es oss 20. 5 4.11 3. 06 6. 59 72.6 2. 49 . 44 2. 14 7.2 
Qo NGO: 40 2 2 dh Fy Sek 20. 5 6. 41 2, 30 7. 09 78.0 1. 99 . 54 2. 62 8.8 
Oma No 23 a. 2). i. ee 20.3 7. 36 . 94 Ooo 58. 7 athe! . 08 .39 1.3 
NS OON Cli ao tp Bee SOc a SLR 19. 0 7, 82 47 Gt 82. 7 1. 38 . 27 1,31 .o 
ORVOPNOf G52. 2 te. 5. EE 19.8 13. 49 5. 67 8, 45 91.8 ike) 720 1; 12 5. 0 
Ohio NO. 77-=— + te. Se 20. 3 15, 88 2. 54 8. 68 95. 2 . 44 won 1.55 7.0 
ORIGIN O88 Se Ao ee 20. 4 17. 54 1.51 8.77 96. 3 5B35: . 58 2. 82 12.7 
ONTOMN On Oieeeh oo ees ae 20.4 | 19,02 1. 48 Soe 95. 9 . 38 . 56 2.72 12,2 
INGEN O 104 ees =e ed 20. 3 20. 36 1.34 8. 32. 91.3 . 80 03 2. 58 11.6 
ORIGIN GO. S49 so oo ote Se 22.3 36. 12 15. 76 7. 74 88. 0 1. 04 . 30 1. 46 20. 0 
rkO) CO 5x6 henna agree caps nee eget eB 21.1 36, 34 Pies 8. 29 92.3 . 68 . 79 3. 84 3.1 
ONTVOING AGT 5 beg cee 22.1 39, 27 2. 93 7. 64 86.8 TL, . 44 2,14 36. 6 
5G) a3) em oe ea aint et ae 22..5 389. 65 . 38 7. 43 84.9 1. 32 ar (1s 3. 00 wraee ets 
OMNOMN G44 Opp ee ee oe a 22.6 40. 11 . 46 5. 87 67.1 2. 87 . 89 4, 32 75.5 
WrorGN GO 482 cs 23 ee 22.5 | 40.72 . 61 6. 91 79. 0 1. 84 . 74 3. 59 62.8 
SO TEROMING, AS et eee Se take 22538 40. 99 ib 6. 76 77.0 2. 02 262 3.01 52. 7 
TONG Hae eoss eee 22. 0 41. 14 ~ LO 6. 82 Ulead 1,01 . 64 Saal 56.9 
Oo NG pues ston eee eee 22.9 48. 86 (Ga?) 8. 43 97.1 . 26 2 Oe 2.02 52.6 
OimorNOAGLee a} .o. oe 22.9 50. 14 1, 28 8. 42 97.0 BIA . 56 ye? 56. 8 
Ole Noro e 22 sas teateo 22.6 50. 48 . 34 8. 21 93.9 Roo . 69 3. 35 70. 0 
ORONO. S042 bon oe eee oe 23,3 65. 73 15. 25 7.31 84. 7 1 aBY) +20 1.21 45.1 
Cirmiberiand 225 o.sees eck 24. 5 66. 11 . 38 6.77 80. 1 1. 68 Roy? 1. 55 10. 6 
OHIOHINO 92020) 22 et 23. 7 66. 48 sot 7. 46 87.0 Peeell 438 2. 09 92.5 
IRONMESGOO = 2 06 A eee ot boos 26/01) 66.53 . 05 7. 58 89.8 . 85 . 16 . 78 12. 4 
OMIGWNGN 958. 22 sos sees e ee 24, 1 67. 11 . 58 7.18 84.3 1. 33 . 30 1. 46 87.6 
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Taste No. 1.—Average results, by months, etc.—Continued 
OCTOBER 1-15, 1914 





Mean time of 



































ae : Total oxygen 
flow, days, Initial dissolv xygen 
a ed oxyge Loss demand 
during 
beat hours’ Quan- 
é i tern- incu- tity 
Sampling station pera- se - é Satu- | bation ' nee 
a- er | ration at parts 
shee Pitts- | tion — fete cent | deficit | 20° c. | Parts per mil- 
burgh | next a Hon satu- | (parts | (parts ilion lion X 
above ration [per mil-per mil-|"™""°"| thou- 
lion) | lion) . sand 
second- 
feet) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 
AlleghenyaNno. 7.9525 --i-8 ATO Neo 5S ee a oe eee 8. 63 88. 7 eat 0. 23 142 0.9 
Monongahela No. 12_-------- LOM = 4 A 5. 97 64. 8 3, 25 . 86 1. 70 1.4 
@hio No: 362 =. 2222-22288 19.0 2a 2, i ota 42 8. 63 . 20 | 14d 252 
Ohio Not Mee. se et te 18.3 8. 34 6. 23 6.37°| 67.2 ee il ae if 1. 3t 2.2 
Ohio Norigee. She. ete 18,8 12 12:78 4, 37 6. 79 12.3 2. 60 .19 . 92 1.5 
Obio No, 2322). 8). 85 < -saee 17.9 | 14.60 1, 89 6. 73 70. 4 2. 83 mali . a8 9 
Boavet. 2. lL ease ae ' 16.0 | 15. 54 94 7. 94 79. 8 1. 62 . 42 2, 04 uf 
@iiio Nio; Gist: ss: Stee c. -tee 17.7 | 25.95 | 10. 41 9. 00. 93. 8 . 50 . 44 2. 14 658 
@iito" No: Wiss. eee ae 18.2 | 30.69 5, 02 | 8. 70 91.6 . 80 Bae 1. 60 4.0 
Oho INOVSS.2 ca. te 55 Ee 18.6 | 33.70 2.78 8. 45 89. 5 . 98 34 1. 65 4.1 
QbIO NO: B7e.- Ba: ee 18.2 | 36.36 2. 66 8. 48 89. 2 1.02 . 38 1, 85 4.6 
Obio No; 4042 . 2258-2 ee 18.1 | 38.74 2. 38 7. 83 82.3 1. 69 . 29 1,41 BsO 
Ohio No. 34052. 3. 55.- ---2 RE 20.2 | 68.38} 29.64 8.17 89.5 96 . 38 1.85 8.2 
SClOlO.: 2 bea eee ee 19.4 | 68.83 45 8. 68 93. 5 | 60 78 3. 79 1.9 
Obip No, dle: 22@.2- ee 20.1 }« 78. 25 4.42 8. 88 97.2 27 Ey 1. 80 13.0 
Die@king . Tee 2. Po: eae 19.3 | 74.12 87 8. 04 86. 5 1. 26 .7 3. 84 av. 
Ohio NO. 4758.42. 22-3 20.:6.45. 75: OF 95 4.75 52.5 4,31 132 6. 41 54.0 
Ohio: No. 4828. 24. bt=2..- 2 ae 18.5 | 76; 27 1. 20 6. 82 71232 2. 62 . 84 4. 08 34. 4 
Olio Not 4888- _ 2283..-. =. 20.1 | 76.69 42 6. 82 76. 0 233 . 82 3. 98 33.5 
@lilo-No}-4028 2. 88-32. = ee 20.0 | 76.88 19 7, 03. 76. 7 2. 14 . 80 3. 88 35. 2 
Obi Nol GO8E 22-885... -2 Be 19.6 | 94.68 | 17.80 9.38 |. 101.5 | —.14 <1 2. 48 25.8 
Ohio Norélies] cee. .2.2 38 20.8 | 98.01 3. 33 9: 38,| 103.9 | —.35 “3 3155 37.0 
@btio No: Gis. osae5. 222 Soe 20.8 | 98.71 . 70 9, O1 99.8 02 BiG 3156 37.0 
Qhio No; 904_..2-42-.---.-12 DAs 27, c 29. 00 8. 31 93. 0 63 . 42 2. 04 48.3 
@umberland: ..- i-25-..-.. 22 21.5 | 128. 2 . 48 GW 80. 4 74 .47 2. 28 7.6 
Ohio No, O20: ..1.94...--2e 21.2 | 128.6 42 8. 33 93. 0 63 . 44 2. 14 57.8 
Ponpessee. es: 22d. 2-52 58 21.6 | 128.7 . 06 8. 35 94.0 54 OL 1. 80 Bont 
Ohi No: 9332. 2.48... 2. ae 21.4 | 129.4 «Ft 8.03 |. 89.8 90 31 1.51 59. 4 


MONTH OF OCTOBER, 1914 

















@ini0 No. 4612. . babes. 2 ee WB-O4e 56 1e- c= 9. 00 94. 3 0. 54 0. 66 3. 21 40.7 
Wicking J. 2h Ss css sole ee 16.6 | 57.19 0. 47 8. 88 90. 5 94 1.05 5. 10 18. 4 
OWING. 47be o> sts oo Be 17.8 | 57. 65 - 46 6. 67 69. 6 2. 91 1. 57 7. 62 128.8 
Oo NOG. 4826. .csheces2-- 17.9 | 58.35 . 70 7. 88 82.5 1, 68 1. 06 5. 15 87.0 
Qii6 No. S832 _. 2 .b64552-2 17.4] 58. 64 . 29 7.79 80. 6 1. 67 1. 04 5. 05 85. 3 
QiIOUNG, 492—= 24 es = 5 ESE 1714) 58, 79 15 8. 36 86. 1 1. 36 1, 21 5. 87 106. 3 
Oyo oO, S08e = tse ses -2 See 18.0 | 67. 40 8. 61 8. 73 91.6 oe silk 44 2.14 47.9 
Oitio No. 61lle 2: 22s. se 8 18.1 | 68. 73 1, 33 8. 92 93.8 . 60 . 54 2. 62 58. 6 
GINO UNO. G19L = snes: ast ee 18.0 | 69.08 . 35 8. 64 90. 6 . 90 57 2utd 62. 0 
Bi RRR cn en oC BS SOI LIPS SE NL SN I a A ee le 
MONTH OF NOVEMBER, 1914 
OURO UINGHIORS: oo. 22 = gaa iC oalty i Slog ee alas tee 11. 94 99. 8 0. 02 1, 90 9. 23 105. 2 
ROICRHING J Gina c oes = cicesin ee 7.6 | 37.74 O64) 122 1) 1061 = 2.20 | 10.69 2.5 
QIN GUA Oe. 2 Set So 8.3 | 38. 28 .54 | 10. 69 90. 7 1. 10 2.44 | 11.85 139. 8 
PME ING. (abies. Soc eee 8.3 | 39.11 .83 | 11.41 96.8 . 38 2.09 | 10.15 119. 8 
OMNOUN Oe goes & ou encore 8.6 | 39.43 .32 | 11.03 94, 3 . 67 1, 98 9, 62 113. 5 
ONTO INGSAG ose ee a eke 8.5 | 39.60 .17 | 10.50 89.5 1, 23 2.10 | 10.40 131. 0 
@IMOING: Wee = ers ee 10.0 | 50.04) 10.44] 13.46] 101.1] —.18 1. 47 14 90. 6 
UO ING: Wiles oc dae 2 deen 9.9 | 52.08 2.04] 11.75 | 103.4 | —.39 eo I ee ae eS ee 
CHOWN Gs GLO 2 cece So = + 4a sie 10.2} 52.60 ae}, LL 99. 6 | 12 1. 46 7. 09 87.1 
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Taste No. 1.—Average results, by months, etc.—Continued 





MONTH OF DECEMBER, 1914 





































































































Mean time of 
itt ‘ Total oxygen 
x 
por ag Initial dissolved oxygen _— Homan 
during 
a 24 
Mean hours’ Quan- 
water F ; 
- : tem- | incu- tity 
Sampling station pera- 4 i sa a oe a 
ta- er ration a parts 
ra Pitts- | tion get cent | deficit | 20° C. * ars per mil- 
“>. | burgh | next pe aaa satu- | (parts | (parts milli on lion X 
above ration |per mil-|per mil- thou- 
lion) lion) sand 
second- 
feet) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 
MOO AGL teehee ete Lee 4.7 toh (38) 74 eee eae 11. 07 79. 7 2. 83 1. 56 7. 58 623. 0 
Went tel Pere tenet cree ape Sel od ee 8. 82 0. 13 OR Be keine |e 2. 44 11,85 12, 2 
OHIOONO3470 >. 22. ced egece 4.7 9. 96 14 11. 38 81.9 2. 52 2750 12.14 | 1, 070.0 
GHTOONOUASD ! 25 Foe lee 4.6 9, 13 vere 11. 50 82. 6 2. 43 2. 48 12.05 | 1, 063. 0 
CHIG ING. 488) 22s een te 4.9 9. 25 .12 11, 39 82. 4 2. 44 2.22! 10.78 950. 0 
104 S02 yaa Se cigie deel ae G8 One Ieee 4,2 9. 27 . 02 12. 92 99. 0 14 2. 46 11, 95 25.8 
GIO ING 2492 = ee 4.6 9. 33 . 06 11. 49 82.8 2.44 2051 12.20. },.1, 103; 0 
OhTOING 508 4 Bo 4.0 11. 93 2. 60 11. 78 89. 7 1°35 1.70 8. 25 796. 0 
OHIO NG GUE. he Sey) 12. 29 . 36 11. 84 89.5 1. 40 2. 06 10. 00 965. 0 
Oiig NOs Ole eee 3. 6 12. 44 715 12, 20 91.9 1. 07 PPS 10. 45 | 1, 008. 0 
MONTH OF JANUARY, 1915 
OiMoywNow4ol ts te ee 12 (alee el es ei 12.78 90. 3 AEeoa 1. 46 7.09 | 1, 219. 0 
yh ba gh Ser Ae el ot Ba, ese 6. 52 0. 10 5 Bef) |) Aa pte as es Ft ala a 1, 24 6. 02 55. 8 
OWIORNO Sigs ge ee 1.5 6. 59 . O07 13. 30 94,8 ne 2.12 |. 10.380 | 1, 893.0 
Ohig Nes 482. Se ee 1.5 6. 71 is 13. 26 94. 5 Berar 1. 97 9. 56. | 1, 756. 0 
OlIGANO. ASS 4s = = se Se 1.4 6. 79 .09 | 13.18 93. 8 . 89 1. 84 8.92 | 1, 640. 0 
IVA ra ene ge 1.0 6. 81 302 |. 13:38 93. 8 . 90 2.08 | 10.02} 39.0 
Oi NO. 4920 ee 153 6. 83 . 02 13. i 93. 0 1. 00 1,91 9, 27 | 1, 740.0 
OinowNO. 509 = 2. oo sees aG 8. 46 1, 59 12. 80 91.4 1. 20 1, 89 9.18 | 2, 036.0 
Olioy NOn OI dja ee oo oS 1.5 8. 66 | 6274 12.74 9028 | 2629 1. 99 9.65 | 2, 140.0 
Olio INORGI9_. 2.2 1.6 8. 75 ld aS TS OV 90. 5 lao 2. 02 9. 80 | 2,076.0 
MONTH OF FEBRUARY, 1915 
Olio Nor 4s6l= sich te P| 65.23" 223 Bee 12936 91.0 128 1. 09 5. 30 | 1, 181.0 
igilos Via Pe ee 3.4 6. 24 0. 01 12. 16 91.0 1.18 1.10 5, 34 61.4 
NVC RUT esata el ah Ue ete PR 6.32 . 08 dl i pet aol agen fs ea . 64 3. 11 3a. 
OMoINO 64/622) = Se a 30 6. 41 Sola 12. 78 95. 5 . 60 1. 60 7.76 | 1, 900. 0 
WIM OING 480 8 2 oes ts Sa 3.2 6. 51 Anite) 12. 76 95. 1 . OO 1. 48 7.19 re 760. 0 
LOIN OFA SOS oer Fee ote 3. 2 6. 59 . 08 12. 73 95. 0 . 68 1, 46 F208 \ 17385, 0 
Whe Evils te Coes oes ta 4,1 6. 61 . 02 12,38 |}. 594.5 ci2 1, 58 7. 67 136. 1 
OO NO? 492 ets no eo Bae 6. 64 2138 12. 65 94,3 se 1557 7.62 | 2,004. 0 
Ohio NGs698—. ext a eo 3.0 8. 07 1;.43,| 12,41 92. 2 1. 07 1.61 7.81 | 2,412.0 
OHIO WNOS6 bleh a9. sb pee 3.3 8. 27 . 20 12:51 93. 7 . 87 1. 48 1.19.2, 220. 0 
Ohio NGw619 5 se cd ee 3.3 8. 37 . 10 12. 47 93.3 . 91 1, 44 6.99 | 2, 158.0 
MONTH OF MARCH, 1915 
OIOUN Op 461 eee 2 Laer 3.9 Ga26> | ce lsen ae 11.99 91.1 Ll: 0. 98 4.76 378. 0 
PAGEIOUVEIATN = ero. Sask 4,2 9, 28 0.02 | 12.1% 92. 8 . 95 1. 02 4.95 3.1 
HACHAT Ciuce 2 AAU sontaee oe uk Spee | tk te 9. 39 Beal el pete’ COG ees ee fee pe . 80 3. 88 16. 3 
ONG NOe4 b= 22 sam... 4.5 9. 53 14 | 12.33 92. 2 . 64 1. 48 7.18 604. 0 
OO MO nO 2 ieee eee eels eee 4.5 9. 69 . 16 12, 33 91.9 . 64 1, 29 6. 26 527.0 
Ohio “N Gc488 = -) yee ee 4,5 9. 81 ld) 2329 94, 7 . 68 1, 29 6. 26 527.0 
VETTE eres.” Loeb eters SE Le 5.2 9. 83 . 02 inearal 92. 0 1,03 1, 22 5. 92 17.9 
OTLO UN On 492 ab eer a ea 4.4 9. 89 06 | 12.19 93. 8 . 81 1. 38 6. 70 584, 0 
OHIO NGG 508 <= eee nce ee 4.7 12, 45 2. 56 12;21 94. 7 1, 69 1. 41 6. 85 676. 0 
OMG NOW6lleas564. bee. 4.9] 12.80 38 | 12226 95. 5 1, 57 1,27. 6,17 609. 0 
Ohio NGA619. owes ek gs. 4.9 12. 95 15 12,22 95, 2 . 61 1,15 5, 58 551.0 








ag 


TABLE No. 1L.—Average results, by months, etc.—Continued 
MONTH OF APRIL, 1915 
































Mean time of 
athe Dae Total oxygen 
5 Init s 
AoW. apye ial dissolved oxygen oe demand 
ates WEES 1.7 ane during 
2 
Mean , Quan- 
water he oh tity 
: : tem- atu- p units 
Sampling station pera- Bia. =a, Pee Prstion Babios iia (parts 
ture, | Pitts- | tion |... mi-| cent | deficit | 490q,) per [per mil- 
°C, | burgh | next p Te 4 satu- | (parts ( ets mili on lion X 
above ration |per mil- = ait: thou 
i . lion) Pp lion) sand 
bate second- 
feet) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 
Oto No 46ie 2 beat ee 12.3 12. 46clens oe 10. 00 93: 0 0. 76 0. 58 2.82 112.9 
Little Miamiocbeoscce: locos 14.4 | 12.49 0. 03 9. 92 96. 5 36 . 68 3. 30 .9 
OU Cl hah eee ene ont pen eee see: FEE SE Se 12, 65 . 16 ORR QS yeaa | eee . 60 2. 92 4.0 
OHIO UNG: Aue =e ee ee eee 12, 4 12. 84 .19 10. 19 94.8 . 05 1.18 5, 12 238. 0 
Ohio Norse: 2 25 5. ee aes 12.3 13. 10 . 26 10. 20 94.8 . 56 . 98 4. 76 198.0 
@inior Nia. 488. eee. ee 129 0))|, 1827 Ue: Onl S 93. 7 . 69 . 98 4.76 198. 0 
Vig eeriaietee: Se fe 0 ene ee 14.4 13. 30 . 03 10. 60 103. 1 132 1. 40 6..80 12.9 
Ona NGr4e2 = se ee 119.) 13537 207 -| ~1OxD5 92. 6 81 1.18 5. 73 249. 0 
GBIGENGN O08s ome ee Soe 12.8 17. 58 4.16 10, 41 97.8 | . 24 79 3. 84 185. 3 
OMOUN OOM =e ee 12.8 18. 25 ae? 10. 83 LOT 5 cal Ire 1. 03 5.00 241.5 
Ohio NGr619. 4 ee =--| 12.8 18. 46 a A | 10. 68 100. 3 —~ 03 1. 09 5. 29 255. 5 





TABLE No. 2.—Factors in resultant oxygen formula (1) obtained from data in 
Table No. 1 (by months) 


Notation: 

Da=Observed dissolved oxygen content of the river at the upper station, in terms of parts per million 
of deficit below the oxygen saturation value at the mean river water temperature between the two stations. 
Ds= Observed dissolved oxygen content of the river at the lower station, in the same terms as (Da). 

La= Corrected initial oxygen demand at the upper station, in terms of parts per million of oxygen. 

t= Mean time of flow from the upper to the lower station, in days. : 

T= Observed mean temperature of the river water, in degrees Centigrade. 

K:1=Coefficient of deoxygenation at the observed mean river water temperature. (Calculated from a 
value, 0.1 at 20° C.) 

Ke= Coefficient of reaeration, calculated from the other terms. 
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At temperature of river water 


























Mean 
Se ee aE oe ee ee velocity 
Month by ies ee ay La ao e 5 2 
(parts | (parts | (parts t Ay K K é oe 
per we - “per | (days) | (° C.) 1 ; weeps 
million) |million) |million) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 
Station 8 to station 11: 
IMA VETO 4 40d 0. 89 0. 63 2. 54 0. 21 15-9 0. 083 0. 990 1, 38 2. 48 
SuNete =. heh ee: 3.13 1. 09 1.91 1. 07 24. 5 123 . 544 35 . 
JUL VRE. Pe oh ee 4.00 1. 32 23 1. 81 24.6 124 . 380 25 29 
ANISUStao Se cso et 5. 53 1.95 4, 22 2.73 24,8 125 . 276 17 19 
September_s.3--52-0 ee 6. 28 2. 49 3. 59 3. 06 20.8 104 197 18 18 
OCchwIaloL ate ee 8. 63 Seid 3. 30 6. 23 18. 7 094 105° 12 09 
Station 11 to station 19: 
ING HG Tee) 1 GR ee eee oe . 63 a 4.89 .18 15.9 083 245 34 2. 50 
Witenes. Oe ee ee 1, 09 1.14 22:20 Bilis 24.5 123 192 12 59 
July a ee 1, 32 1. 04 2. 54 1.36 24.1 121 289 20 Be 
AtIOUSh: 2:20 4 2 1.95 1. 60 3. 44 1.98 24. 4 122 217 14 ‘22 
Sapvusmipen. Oe. gases 2. 49 1.09 3.14 2, 30 20. 5 102 144 14 19 
Octwials. 2 2e.t 23 a ee Ps apt 2. 60 1. 74 4. 37 18. 6 092 052 06 - 10 
Station 23 to station 66: 
Misys Old: 38... Bo 85 1. 22 5.16 . 62 16.3 084 200 30 4, 20 
dametwe. See ot ete 1. 33 . 93 2. 94 1. 90 22.9 115 330 25 1. 36 
Vulivese... Seek. Fe 1, 24 . 94 2. 65 3. 00 24.8 124 230 15 86 
PAP oe Rh coisa etebs 1, 88 wt5 3. 69 4, 90 24. 4 121 260 ive: > 53 
September. ees 3.77 . 76 2. 56 6. 14 20. 0 100 190 19 . 43 
Oeil by SRE. 2. I. 2. 8&3 . 50 15. 74 LIVSO 18.0 091 353 43 | 
Station 66 to station 77: 
IVI EVa Ola ee ec ee we 1. 22 12 4.89 .18 15.9 083 245 ung 3. 64 
Ut Osere ae et . 93 . 46 3. 08 1. 02 22 116 801 59 7 
VtDiver ee eae see ee . 94 . 64 2. 14 Lis 24.8 125 436 28 44 
AUIRUStaee ee cae Sal's 75 . 49 1. 26 2. 47 24.4 OB e268 st yi .31 
Soptembperseesecces.ce.ce . 76 . 44 1. 92 2. 54 20. 0 . 100 . 334 . 33 231 
OGUEIGee ee ere U ee se . 50 . 80 3. 34 5. 02 18. 0 . O91 . 189 APS} .15 


| 
| 


3” 


Taste No. 2.—Factors in resultant oxygen formula, etc.—Continued 





At temperature of river water 
































1S Mean 
Ka (at ye gre 
Da Ds La 9 (at | of flow 
Month (parts (parts (parts t T K K 20° C »)) (feet 
per per per | (days) | (° C.) t 3 per 
million) |million)|million) second) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 
Station 77 to station 88: 
Mayen 1914.4 Akos acl 0. 98 0. 89 6. 30 0.18 16.8 | 0.086 | 0.802 07 3. 68 
Fiimere 8} ee oe ae . 46 . 389 3. 55 . 65 23.7 . 118 1. 005 71 1. 01 
POON. ober 8 Fee 5 54 -b2 2. 56 1. 01 25.0 . 126 . 529 . 33 . 62 
SUMS tort es eel ot .49 S/ 2. 15 Tbk 24,5 . 123 . 245 16 40 
September: . 2... ek . 44 235 2. 87 1.51 20. 4 . 102 . 716 .70 40 
IS] jal aad 9 iy ee eon A . 80 . 98 2. 29 2. 73 18. 4 . 093 . 144 . 16 wae 
Station 104 to station 349: 
Nia, thee 1. 50 1.31 3, 80 4.03 16. 7 . 086 . 180 ! . 24 3. 2 
UG: cet 3 EA A A . 96 . 88 19.1 10. 09 24. 4 . 124 . 268 .18 1, 48 
EN oe Ca i ea a . 92 1. 28 41.6 12. 79 25.8 Filo . 202 .12 1.18 
LSE TE Ra Sa Se a . 95 1.18 | 115.8 16. 62 24.8 +225 mee) .14 91 
September__.--..--.--.21 . 80 1. 04 dae. 15. 76 21.3 . 105 .172 15 . 96 
GTS. 2S re 1.69 .96 | 618.2 29. 64 19. 2 . 097 . 200 . 22 . 51 
Station 349 to station 461: 
1 ENTS G 1 he Eo eee ae Pst 1, 42 2. 66 1. 76 bid . 090 . 128 16 3. 90 
Silde-s Ste a ee . 88 21 3. 90 2. 92 25.8 . 180 1. 140 . 65 2.36 
CEU 3 Se eee Pe 128 35 2, 44 2. 96 26.9 . 142 . 530 . 28 2. 32 
MOsust sl! tiei he Firte 1.18 . 89 orl 3. 46 26. 0 . 137 . 258 15 1. 98 
September..-_2.........5 1. 04 TZ 2. 84 3.15 2252 rade a LOg 18 ‘2.18 
OGEI-19 5 SES es . 96 64 4.08 4, 87 20. 2 . 102 . 340 .33 1.41 
Station 475 to station 482: 
Way; 1914. £0230: gi 1. 54 1/55 4,02 543 77 1 . 088 5229 . 29 8. 35 
2 ae ES ck" hae HER 1. 84 1.46| 4.77 . 48 26. 0 5132 J 2. 820 “a2 93 
Fulby fac2 2} i nee cal ES 2.02 1. 34 4,21 = +) 66 27.0 . 189 . 650 . 33 79 
SARURP TIS Gee). ce A 3220 1. 86 4, 20 . 59 26. 0 2 i322 . 460 . 26 . 65 
September_...._-_______. 2. 87 1, 84 4, 28 . 61 22.6 . 113 . 514 . 40 3 
OGEO DOG: 2) ae el 2.91 1. 68 6. 75 . 70 17.8 . 090 . 601 te . 63 
November. if sieye- 1.10 .388:| 11.60 . 83 8.3 .058 | 1.750 3. 98 . 54 
Hecem bere. = 2 2.52 2. 43 12, 24 Ae, 4.6 . 049 . 345 . 96 2. 69 
Joniary, 1915... wie) v3 rig 10. 00 12 125 . 043 . 365 1,14 3. 83 
IPGDEUAIY sa = Bee . 60 . 65 7. 50 . 10 3.2 . 046 . 202 Od 4, 29 
Moaypelna te) fd) -ziesh . 64 . 64 6. 76 . 16 4.5 . 049 016 1, 43 2. 74 
LG SE SE See £43) . 56 5. 35 . 26 12. 4 . O71 . 612 1.09 1. 73 
Station 482 to station 488: 
TG o3\ A) 2 See ees 1.65 1. 60 Sol .10 17.0 . 087 . 040 . 05 8. 61 
Jane sath otace sors7 1. 46 1. 83 4.08 25 2527 VASO ala. Seeseb ies Pe ey 1.47 
| (1 Se a ee, ee 1. 34 1. 42 3. 76 26 27.0 1388 . 258 13 1. 40 
maIsaStey eG Pe A. 1. 86 2. 00 3. 47 29 25: 1 126 . 118 07 1, 23 
DEMUGHI DOR o.= a = eS 1. 84 2 O2 SOL. 27 22. 4 2 . 060 05 1. 34 
October tet. x cities 1. 68 1. 67 5. 20 29 17.6 090 . 285 35 1; 23 
INOVOMIDOE 2-0 . 88 .67 | 10.10 32 8.4 059 . 858 81 1,12 
Pecenther ee! sey evi sb 2. 43 2. 44 11, 63 12 4.8 050 EOTT 39 3. 00 
January,<al9tb. 2.2.2. TE . 89 9. 35 - 09 eee 043 . 033 10 4,11 
Pebritary. 12580232202 0u . 65 . 68 7.138 . 08 3.2 046 . 265 78 4, 57 
A aR aaa Maley . 64 .68 | 6.29 .12 4.5 049 | . 265 71 3. 00 
Aprils 4233 Sos . 56 . 69 4, 87 1% 12.2 070 . 075 14 2.12 
Station 492 to station 598: 
Weary, (iQ. ste" sree ais 1. 67 1. 61. 4, 74. 1. 98 17.2 . 089 2213 27 3. 26 
TED a Ere ee See 1. 90 82 15. 38 6. 06 Ii a 126 0382 32 1. 07 
Jolytes th errr. 7 ek 37 24, 80 7. 26 26. 8 135 1. 090 57 . 89 
NGG Gt FF SR DE 2. 15 41 27. 30 8. 67 25-5 128 788 46 74 
September sin2 2252 2. 2. ry 1 a0 . 26 11. 06 Mal? 22, 4 112 . 761 61 84 
Wei Debora e eee 1. 36 .8l 9.18 8. 61 17.6 090 . 247 31 ihe 
Naperiberhi 704 ware oenter: Per ee Pert tees gg Mes ol eer yee a rigiae: Pepe 
Weeemher.. = ..-.225, 525 2. 44 Too 11. 60 2. 60 4.3 049 . 380 1. 05 2. 48 
Janwary, Wirt zccccehs. ; 1.00 1. 20 10. 80 | 1. 59 14 043 . 807 99 4, 06 
TYE) 9) 6k FS 1 aay io) 1. 07 8. 41 1. 43 ied 047 . 200 59 4. 61 
M arch ee eer a 2 £8 1. 38 1.41 7.98 2. 56 4.6 050 . 234 65 12.H8 
April 2 2 a ee er ee as . 81 . 24 6. 69 4.16 12. 4 071 1. 050 1. 88 1. 56 
Station 598 to station 611: 
IMWay 1914. -* 7 1. 61 1. 64 3. 39 . 28 17.8 090 170 21 2. 94 
dUNe@x2i cho one At . 82 . 86 4,97 .95 25.6 129 613 36 85 
Wy atec ae ROR ond es ot . 59 5. 41 1, 20 Zlal 139 925 47 68 
INMUIOUISh bas ent FS 41 . 28 2422 1, 44 26. 5 135 792 43 56 
September... ks . 26 ae. Se Ld 1. 28 22. 9 114 | 1.140 87 63 
Octoberes 2% cies. eeu ss3 . 81 60 4. 68 3783 18.0 091 446 54 60 
INOW OIRVEL.<.. .-2200--5 —.13 —.39 8. 90 2. 04 10. 0 AOGOe saan soe anh ommlag brome ee 
December _ -.2:2..2.----- 1.35 1. 40 9, 20 . 36 3.8 047 264 15 2. 24 
MAMA, POLO. ooo ce ae 1. 20 1, 29 9. 57 21 1.6 043 185 58 3.73 
Pebrivary2t clvwceteerl: _ 2 1. 07 . 87 7. 58 .19 Bae . 046 830 2. 44 4.11 
arch RE el tS OM 2S ES asad 1, 69 i eeay 6. 94 . 00 4.8 . 050 295 82 2. 28 
ADF . osc. po easeaesap- . 24 —.18 4,75 $72 12.8 BOTA MC eOsrt Pa eet i cup ee 





1 Results abnormal; omitted. 
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TaBLE No. 3.—Derived values of (Kz) at 20° C. (from Table No. 2) 









































Sep- : No- | De- | Jan- 
Stretch May) June July Au- tem- Oct. | Octo- yem-; cem- | uary Feb- March; April 
(1914) gust ier 1-15 | ber per. [iber | (1915) ruary 

Sates 5 a Lt. 381! 0.338 710 25. | O17 1 Oe TSO EDM RA Re La AS Daa Se ae SS ek oles eee 
hE) {Sa ee eer tt 34 Az ReaD) 14 14 (0 Gielen Dt, omen Rohl fe ERS | agen alee 
23-6b..2 ene ee ee 30 . 20 Pi 9) ay, 19 7 al ae es ae, Sire | 7 ee | SP ae | genes SS & fe 3S 
G5-Teih. 1. GD 23 111), CBO). 28 17 {23.33 Sih (i oe ey ee eet ee in 
Ufa eee oa Gee este 1. 07 ire 38) 16 . 70 GS es ee eee call Seeeeeah tes 44 metal on meen alee eee 
MOGS SAO Ae ar oe 24 . 18 Fp 14 Py t's) DO se A ee ) ee Oe ol eephe = iee e 
OA OAGN = tee 16 . 65 . 28 15 we p23 | a Wee Oe ee a ee il a Pe RE MM 
A748 2. 4. SM ot. LES 29 exe P'4 £33 26 40 Wee: 0.73 | 3.98 | 0.96 |-1.14-}. 0. 61 1. 43 1. 09 
ASDA48S. SANS. sae 05} |Sae . 13 07 06 eee 35 81 39 10 78 ve 14 
PA TS me es ey Ee 27 32 57 . 46 OR pe Se Ly i eee ee 1. 05 99 59 65 1. 88 
Ea; hoc Op UE Ska pee 21 36 47 . 43 FeV jie eee oe a eee 76 58 | 2. 44 BO aus 




















THE OXYGEN DEMAND OF THE OHIO RIVER 


The results of oxygen demand observations in the Ohio River and 
at the mouths of certain major tributaries, which are given in Table 
No. 1, columns 8, 9, and 10, are presented in terms of the 24-hour 
loss of oxygen on incubation at 20° C. (column 8), and also in terms 
of the total oxygen demand as calculated from the 24-hour figure 
by taking the value of the deoxygenation coefficient (K,) as being 
0.1.(see p. 7), whence the 24-hour demand represents 20.6 per cent 
of the total. ‘It may be noted in this connection that this value of 
(K,) defines a deoxygenation curve coinciding exactly with the 
“relative stability” curve given in Standard Methods of Water 
Analysis of the American Public Health Association (1920 edition, 
Table No. 15, p. 70). In column 10 of Table No. 1 the total oxygen 
demand figures of column 9, in terms of milligrams per liter, have 
been weighted in each case By the discharge of the river, in thou- 
sand sauatidls feet, observed at the particular sampling station co- 
incidently with the laboratory determination. The purpose of these 
weighted figures is to show the total amounts of unoxidized organic 
matter carried by the river at various points, corrected to a common 
basis of dilution; in other words, to eliminate changes in dilution 
as factors in masking real changes in the oxygen demand of the 
river between successive stations. Strictly speaking, the unit chosen 
for this purpose is one of rate; that is, it measures the quantity of 
biologically oxidizable matter (in terms of oxygen demand) carried 
by the river past a given point in each unit of time. Thus, the fig- 
ures in column 10 may be reduced to terms of grams per second 
by multiplying them by the factor 28.3. For convenience, however, 
these values will be designated as “quantity units.” 

Corrected in this manner, an increase in the number of “ quantity 
units ” of oxygen demand as observed between two stations indicates 
that added amounts of oxidizable matter have been brought into the 
river at intermediate points either in the form of sewage and other 
wastes discharged directly into the stream, or through the medium 
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of tributaries and local surface drainage. A reduction, on the other 
hand, indicates that progressive satisfaction of oxygen demand has 
more than counterbalanced the effects of intermediate inflow in the. 
river stretch in question. Ifthe oxygen demand data be interpreted 
from this standpoint, some highly significant facts are disclosed as 
to the sources and behavior of oxidizable matter in the Ohio River. 

In certain stretches of the river, relatively undisturbed by inflow- 
ing pollution, marked and consistent reductions in its oxygen demand 
are observable during periods of settled flow conditions. Among the 
shorter stretches, this tendency is particularly noticeable in those 
extending from stations 3 to 11, below Pittsburgh, and stations 475 
to 488, below Cincinnati, where decreases in the oxygen demand of 
the river were observed consistently he sees the summer and autumn 
months of 1914. 

Of the longer stretches, which are in general subject to relatively 
large additions of inflowing drainage water, even during low water, 
the only one in which a reduction in the number of “ quantity units” 
of oxygen demand was observed with any degree of consistency was 
in the stretch, approximately 300 miles long, extending from Station 
619, below Louisville, to Station 904, a comparatively short distance 
teow the mouth of the Ohio River. During the period of June 1 
to October 15, 1914, the average number of “quantity units” of 
oxygen demand observed at Station 619 was 81.3, or 2,300 grams per 
second. At Station 904 the average number during the same period 
was 50.6, or 1,482 grams per second. In spite of the large volumes 
of pollution brought into the Ohio between the two points by the 
sewage of Evansville, Henderson, and other cities, and by the Salt, 
Wabash, and numerous smaller rivers, the observed reduction in 
oxygen demand in this stretch of the river was about 38 per cent. and 
it would undoubtedly have been much greater in the absence of the 
intermediate sources of pollution noted. 

In longer stretches of the river farther upstream, as, for example, 
between Station 104, below Wheeling, and Station 349, immediately 
above the Scioto River, the influx of a number of highly polluted 
tributaries, adding a large volume of flow to the Ohio, has a constant 
tendency to mask the true effects of organic oxidation going on in 
the stream, as is evidenced by the increase in the “ quantity units” of 
oxygen demand between these two stations even during the summer 
low-water period of 1914. A similar tendency, though not so con- 
sistently marked, was observable in the long stretch of the river ex- 
tending from Station 492, below the mouth of the Miami River, to 
Station 598, above Louisville. 

In one stretch of the Ohio River receiving almost solely rural 
drainage, namely between Station 358, below the Scioto River, and 
Station 461, above Cincinnati, the increase in oxygen-demand values 
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observed during the study was so marked as to indicate that a very 
considerable proportion of the total oxidizable organic content of 
the river originates in surface drainage from sparsely inhabited and 
unsewered areas. This conclusion is borne out by the fact that the 
Ohio River carried, on the average, a very considerable excess of 
oxygen demand over that which would be accounted for by the total 
sewage contribution of the watershed. At Station 461, above Cin- 
cinnati, for example, the oxygen demand of the river during the year 
extending from May, 1914, to April, 1915, inclusive, averaged 352 
erams daily per capita of sewered population in the entire water- 
shed above this point. From observations immediately above and 
below Cincinnati during the low-water period of June 1 to October 
15, 1914, the average amount of oxygen demand contributed to the 
river by the Cincinnati metropolitan district during this period was 
equivalent to 140 grams per capita daily, which, it may be noted, does 
not greatly exceed the figures for normal sewage, 112 grams and 
100 grams per capita, given on pages 2-3. Taking these values as 
representing the normal contribution of oxygen demand in the form 
of sewage, the large excess of this constituent carried by the Ohio 
River, over that which would be accounted for in terms of sewage, is 
evident. 

It may therefore be said that while evidences of the forces of pro- 
gressive oxidation at work in the Ohio River are discernible in 
certain stretches of the stream and under conditions wherein sta- 
bility of flow exists, these evidences are for the most part masked 
by disturbances resulting from inflowing pollution originating both 
in sewage and in surface drainage water. 

Bearing in mind these factors, it is noW proposed to examine the 
data with reference to the validity of an assumption made in de- 
riving the formulas on pages 18-19; namely, that the rate at which 
oxygen demand is satisfied in the stream is governed by the same 
law that controls its rate of satisfaction under conditions of the 
laboratory. This law has been stated on page 5 to be as follows: 

The rate of biochemical oxidation of organic matter is proportional to the 


remaining concentration of unoxidized substance, measured in terms of Oxi- 
dizability. 


Or, in mathematical terms (p. 5): 


logy = Kt 


Owing to limitations above stated, the oxygen demand values ob- 
served at the various Ohio River stations must be sifted carefully for 
evidences of disturbance from surface inflow before being subjected 
to so severe a test as the one proposed. It is desirable, moreover, to 
have available for the purpose data covering a fairly long period, a 
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year at least. The only data from the Ohio River studies fulfilling 
these specifications are those which were obtained in the river stretch 
below Cincinnati extending from Station 475 to Station 488, a dis- 
tance of, roughly, 13 miles. About midway between these two points 
was Station 482, but since it was located but a short distance below 
Government Dam No. 37, over which the water has a free fall during 
periods of low water, the oxygen demand results obtained at this 
station were to some extent influenced by the presence of small 
amounts of entrained air in the samples and were therefore not as 
reliable as those obtained at Stations 475 and 488. Tor these two 
stations results were available over the full period of May, 1914, to 
April, 1915, inclusive, as given in table No. 1, and in addition, for 
the months of May and June, 1915, during which period certain 
special observations were continued in the river below Cincinnati. 

Considering these monthly average figures at Stations 475 and 
488, 16 was obvious that if the rate at which the oxygen demand of 
the river was satisfied in this stretch were a logarithmic function of 
time, in accordance with the formula above stated, a close linear 
correlation should exist between the observed times of flow between 
the two stations, corresponding to (t) in the formula, and the 
logarithms of the ratio: 


Oxygen demand at Station 475 
Oxygen demand at Station 488 


corresponding to (log =) in the formula. Employing the Galton- 


Pearson coefficient (r) as an index of this correlation, a value of 
() =0.85 0.043 was obtained from the observed results. Since a 
value equal to unity would denote perfect correlation, it is apparent 
that the value obtained indicates a high degree of correlation between 
the two variables as observed in the river; so high, in fact, as to leave 
little doubt concerning the agreement between actuality and theory 
in this case. 

A more severe test of this correlation was provided by an analysis 
of the same data from another viewpoint. Referring to Table No. 
4. in which the results of this study are tabulated, the observed 
monthly reductions in oxygen demand between the two stations 
were converted to terms of percentage of the initial value (column 
5) and reduced to a common basis of loss in 24 hours at 20° C. 
(column 8), using for this correction the formule and temperature 
factors described on pages 7-8. Referring to column 8 of the 
table, it is noted that during the period June—-November, 1914, when 
flow conditions in the river were stabilized, the decrease averaged 
23.5 per cent in 24 hours. At the mean temperature observed at 
Station 475 during this period (21.3° C.) the value of K,, as given 
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in Figure No. 8, would be 0.116, and the 24-hour decrease, assuming 
a logarithmic rate, would theoretically be 22 per cent. For other 
months than the June-November period in 1914, a wide variation in 
the 24-hour figures was found, but it will be noted that during these 
months flow conditions in the river were unsettled, and observed 
times of flow between the two stations were so low that large errors 
of extrapolation would be involved in converting the results to a 
24-hour basis. 


‘TABLE No. 4.—Reduction in oxygen demand between Ohio River Stations AVS 
and 488, as observed and as corrected to a ere basis of time equivalent 
to 24 hours and temperature 20° C. 


[Arranged in order of magnitude in times of flow] 



































Oxygen demand, D bat Ti ff d 
parts per million ecrease between ime of flow, days, Per cent 
eR stations Stations 475 to 488 decrease, 
Month noe hee BOS, 
20° C. (at 
4 At river | Equiva- | observed 
Station | Station | Parts per| Per cent f 
475 488 |’ million | of initial | *™Pera-| lent at | rates) 
(1) (2) (3) (4) (5) (6) (7) (8) 
I. High river stages, shorter | 
times of flow: : 
Mebrilary, [ines race ee 7. 76 7. 08 0. 68 8.8 0. 18 0. 11 ta O 
January, [9152s 26" ee. 10. 30 8. 92 1. 88 13. 4 Zi. . 09 134. 8 
WAY, Lolact nets oe ete 3. 84 Ayu alk wigs 19.0 | . 23 noe 81.5 
Maren. 1005: 247 tase oo ee 7.18 6, 26 . 92 11.4 28 .14 re et 
POCEIMDER, U14s 222 eee 12. 14 10.7 1. 36 2 . 29 . 15 65. 8 
Vie. LOUS2 Mie sao Bae ale 3. 45 2. 86 . 59 17.1 31 . 34 47.1 
INTO; - 1016. 26 Se eeeee es 5. 59 4. 42, lV, 30. 4 exh . 30 80. 4 
April; 1915e.4 +: wa ceecersicr 245. 72 (4.76 . 96 “~t6: 8 43 . 30 50.8 
II. Low river stages, longer | . 
times of flow: 
June; 19904. ge fees FLY 4, 37 3. 60 Ta 17.6 73 94 18.6 
RD Ce a ESE ie 3. 94 3. 30 . 64-] s 16.3 82 1. J2 14.7 
September, 1914.__._-____.. 4, 82 3. O1 1:32 30.3 88 98 30. 9 
STIRS, LON Ae sr tee oe 2d 3. 84 3. 15 . 69 18. 3 98 1, 28 14,8 
October; 1914 -rs_ sient ts 7. 62 5. 05 23 57 33¢7 99 99 33.8 
November, 19142... 2 2. 11. 85 9. 62 223 18.8 115 70 25.8 
Oct 15, 19143 lara 6. 41 3. 98 2. 43 37.9 1. 62 1. 56 25.8 
Mean, June-November, 
i 0 0 ‘aa A 2S | | eer Semeese es |e ee ON SE SS See 23, 5 





Further and even more striking evidence of the closeness with 
which the rate at which oxygen demand is: satisfied in streams fol- 
lows the logarithmic law has been afforded by a preliminary anal- 
ysis of extensive data obtained by the United States Public Health 
Service from studies of the Illinois River during the years 1921 and 
1922. Observations in. a stretch of this river about 70 miles long, 
extending from La Salle to a point a short distance above Peoria, 
where conditions affecting natural purification are in a very fair 
state of equilibrium, were especially illuminating on this point. 
Taking the mean oxygen demand values observed at La Salle (Sta- 
tion 227) as the 100 per cent point, and computing the oxygen de- 
mand amounts observed at three other stations downstream as per- 
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centages of the La Salle figure, these results have been compared 
for the period of December, 1921, to April, 1922, with corresponding 
percentages calculated on the assumption that a logarithmic rate of 
decrease in the demand held throughout this river stretch. The re- 
sults obtained are as follows: 


5-day oxygen demand, 
per cent of Station 











227 value 
; Time of 
Station 

flow Assuming 

ogarithmic 
Observed mateot 
eae 
Hours 

127 (La Salle) qc; 23 hue ee sort tld. dep pete s. -ob ea eins 0 100 100 
RCRR ESE COTA TEN; ieee een oe SW oe Me ogee sR ac ieee 45 84 71 
£79 1 (@ Di eOtne) sake fe aoe 2 Ree aE as eae 103 43 47 
SOC AC ET WAVING) is skye ear Fd it Cals sgn we oa eT, ol Oats Se ae 144 aD 135 





1 Assumed the same as the observed figure in order to calculate values at Stations 196 and 179. 


The closeness of agreement in the observed and calculated values 
at the two intermediate stations, 196 and 179, is the more strik- 
ing in view of the relatively high river stages and low tempera- 
tures prevailing throughout the period of the observations. A 
similar comparison, based on observations during the summer 
of 1921 and 1922, was less satisfactory in view of an observed 
increase in oxygen demand between Stations 179 and 166 during 
this period, probably as a result of the intensive growth and de- 
composition of microscopic plants in this stretch of the river, which 
in the summer presents conditions approaching more closely those 
of a shallow lake than those of a flowing stream. Based on the 
summer observations between Stations 227, 196, and 179, however, 
the observed oxygen demand at the intermediate station, 196, was 
85 per cent of the Station 227 figure, as compared with a calculated 
value of 81 per cent, assuming a. logarithmic rate of decrease. 
Inasmuch as the total decrease in demand between Stations 227 
and 179 amounted to 34 per cent of the Station 227 value, the 
close agreement of the calculated with the observed demand at 
Station 196 is significant. 

In considering the foregoing data it should be borne in mind 
that observations of this kind made in natural streams are sub- 
ject to. inherent errors such as would cause inevitable variations 
in oxygen-demand_ results obtained even under stable conditions 
of flow. The use of such results for the purpose of verifying their 
conformance to a fixed mathematical law, as has been done in the 
present case, constitutes an extremely severe test, to be inter- 
preted rather broadly and with due allowance for the inevitable 
“spread” of the data, Viewed in this light, the evidence presented 
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above affords ample justification for the assumption made in deriv- 
ing formula (1). 3 


THE OXYGEN RESOURCES OF THE OHIO RIVER 


It has already been shown that the total oxygen resources of a 
stream at a given point are, first, the dissolved oxygen actually 
contained in the water at that point and, second, the potential re- 
sources available through the capacity of ihe stream for reaeration. 
The data relative to these two sources of oxygen will be separately 
discussed. 

Dissolved oxygen.—In Table No. 1 the dissolved oxygen values at 
various stations throughout the Ohio and at various times during 
the period of the observations are shown in three ways, viz: 

(a) In parts per million ; 

(6) In percentage of saturation ; and — 

(c) As the saturation deficit in parts per million. 

It is to be noted that the concentration of dissolved oxygen itself 
becomes critically low only during average and extreme low water 
and even then only in comparatively short stretches of the river 
below Pittsburgh and Cincinnati. The conditions for a few miles 
below Pittsburgh approached complete exhaustion of dissolved oxy- 
gen at extreme low water during October, 1914. Indeed, it is doubt- 
ful whether the present comparative freedom from general nuisance 
at this point would be enjoyed were it not for the bactericidal influ- 
ence of the acid conditions in the Pittsburgh district. The lowest 
average oxygen value recorded at station 475, about 7 miles below 
Cincinnati, was 52 per cent of saturation, observed during the first 
half of October, 1914. As far as could be determined, no offensive 
conditions existed in the Ohio River below Cincinnati at that time. 

During low river stages the waters at the mouths of several of 
the tributaries were found to be at a point of lower oxygen con- 
centration than the waters of the main stream. This is notice- 
ably true of the Cumberland and Tennessee Rivers, entering the 
Ohio near its mouth, at which point the main stream usually has 
recovered its HRA Mod oxygen to an extent of from 85 per cent to 
over 95 per cent of its saturation value. 

A study of the saturation deficit figures is instructive in showing 
the rapid rate at which the Ohio recovers its oxygen supply thrdiph 
reaeration, even during low water periods. ‘Thus, as indicated in 
Table No. 1, it appears that at station 461, above Cincinnati, the 
stream is often in better condition as regards reserve oxygen than 
it is above Wheeling, and is sometimes in better condition at station 
904, near its mouth and above the Cumberland and Tennessee 
ee than above Cincinnati. These facts are of special signifi- 
eance in connection with the importance of the reaeration factor. 
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Reaeration—A. direct computation of the amount of reaeration 
in the river, expressed in any convenient concentration unit, may 
be readily made by applying the principle illustrated in Figure No. 
5 and determining the difference between the actual oxygen content 
at the lower station and its calculated value in the absence of reaera- 
tion, as determined from the average oxygen demand for the given 
river stretch. 

The actual amounts of reaeration per mile of river length ob- 
served during various months and in various stretches of the Ohio 
are given in Table No. 5. These amounts are rather widely dis- 
persed, as would be expected in view of the marked variations with 
respect to oxygen saturation deficiency and to physical conditions 
of channel and flow encountered in the several river stretches and 
at different times. The following further analysis of the data was 
undertaken to determine if the reaeration values obtained are con- 
sistent with the known physical features of the various river 
stretches. | 


Taste No. 5.—Average rates of reaeration, in parts per million of dissolved 
oxygen, per mile of river length, observed in designated stretches of the Ohio 
River during specified months of 1914 and 1915. 
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1 Probably in error by a wide margin owing to the excessively long indicated time of flow. 


Because of the varying and complicating influence of these numer- 
ous factors, the reaeration values themselves are only of special 
local significance. A study of derived values of the reaeration coefii- 
cient, however, and of its relations to depth and turbulence, furnishes 
data of fundamental importance, and attention will therefore be 
directed to this phase of the subject. 

The reaeration coeficient-—Values of the reaeration coefficient 
(K,), calculated by the method outlined on pages 20-27, are given 
m Table No. 3. In Table No. 2 are also given the various factors 
employed in the calculation, transcribed from Table No. 1, and in 
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the last column, values of (K,) corrected to a standard temperature, 
20° C., by means of the curve in Fig. No. 7.. The proper symbols for 
the various quantities in formula (1) are given at the heads of their 
respective colunins. 

On referring to Table No. 2, it will be noted that the assigned 

values of (LL) for the two loxgex? river stretches, stations 104-349 
and 492-598, are so high in some cases as to appear. rather incon- 
gruous sie visualized in terms of actual oxygen concentration in 
the stream. This is notably true of the high value of (L,) given 
for the river stretch, stations 104-349, for the period October 1-15, 
1914. On referring to column 5 of the table it will be noted that 
coincident with these excessively high values of (L,) are times of 
flow ranging from slightly less than 10 days to nearly 30 days, 
whereas. the times of flow coinciding with lower values of (Lz), 
both in these two and the other river stretches, are, im general, pe a 
lower order of magnitude. 

In interpreting this apparent abnormality in the (L,) figures, the 
significance of (LL) as defined on page 18 and as further interpreted 
on pages 21-26 should be borne in mind. Being an hypothetical 
value, it represents, in all cases, an assumed quantity sufficient in 
magnitude, if concentrated immediately below Station A, to cause 
a residual oxygen demand at Station B equal to the arithmetical mean 
of (ig) and (LL’,). Where the time factor is long, as in the two 
river stretches above cited, the quantity thus assumed is necessarily 
higher than would be the case if the time were shorter, and it may 
be much higher than the total oxygen demand actually introduced 
into the stream at points intermediate between A and 5, causing the 
divergence of (Ls) from (Li’,). Perhaps the best evidence that these 
excessively high values of (L,) are not abnormal is afforded by the 
general consistency of values of (K,) derived from them, as com- 
pared with corresponding values derived for the same river +tretch 
during months of shortened times of flow and lower values vf (L,). 

Taking as an example the results of (K,) calculations in the sta- 
tions 104-349 stretch, as given in Table No. 2, column 8, comparison 
of the (K,) value, 0.18, obtained for May, 1914, (L.) being 3.8, with 
the value, 0.20, for October, (La) being 618.2, is significant evidence 
on this point. As previously noted, however (see footnote, p. 23), 
errors of observation and extrapolation are most likely to produce 
erroneous results when the time factor is long; hence the results of 
calculations in Table No. 2 for the two river stretches, stations 
104-349 and 492-598, during the months of extreme low water in 
the summer of 1914 are to be regarded as being less reliable, on the 
whole, than corresponding results obtained for the shorter river 
stretches during the same period. These results have been retained 
for further analysis chiefly because of their general concordance 
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with the figures obtained. for the shorter river stretches; likewise 
because they are desirable material for subsequent, correlations. to 
be described later in this text. 

Inspection of the (K,) values, corrected to a uniform temperature 
basis, shows that they are subject to considerable variation in each 
river stretch under the different flow conditions encountered from 
month to month. With a few exceptions, however, their general 
order of magnitude is about the same in different river stretches 
when determined under similar conditions of river stage, indicating 
that the conditions which cause variations in rates of reaeration are 
much more nearly uniform in different sections of the river under 
similar conditions of flow than they are in a given stretch at dif- 
ferent river stages. In other words, those physical factors which 
change with river stage are the more influential upon reaeration 
rates. 

The two stream conditions most affected by variations in gage 
height are depth and velocity of flow, the latter being closely re- 
lated, of course, to turbulence. In Table No. 6 a fairly definite 
correlation is shown between values of (K,) and both depths (H) 
and velocities of flow (V), the depth (H) being expressed as gage 
heights at river gages listed in the footnote of the table. The true 
nature of this relation, however, does not appear until each (K,) 
value is weighted by the square & the depth (H?), whereupon the 
quantity (K,H’) is found in every river stretch to approximate 
closely a simple power function of velocity (V), os the for- 
mula: 


K,H?=cv" (4) 


TABLE No. 6.—Relation between the reaeration coefficient (K:2), the mean depth 
(H) of the stream above extreme low water, and the velocity of flow (V) 
in certain stretches of the Ohio River 


Notation: 

K,.—Reaeration coefficient. 

H =Mean depth in feet of stream above zero points of reference gages, which points 
are taken as representing extreme low water. 

V =Mean velocity of flow of stream in feet per second. 
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(i geese Sagi oi lgag ir go See a SR I NS OS Ore 20 3.6 2. 60 .32 
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TABLE No. 6.—Relation between the reaeration coefficient (K.), the mean depth | 


(H) of the stream above extreme low water,,and the velocity of flow (V) 
in certain stretches of the Ohio River—Continued 
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Station No. 23 to station No. 65: 5 - | c | 
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Mais 1914-0. on ct de dey ee eee. Satay Feats fre. 1. Uy 11.4 144.0) 3. 64 
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INDUS ad se - ce ake 2 oe Se ee eaten ee Mee ec ee a Le 8.2 pal ie. | aii! 
Septem bero.t642 9: - tes, poregi se. cased eee - es } . 33 6.9 VST iE" (3 
OCU ASO oe dee eee ee ets a Se ee Loree 223 | 2.0 .91 15 
Station No. 77 to station No. 88: 
nv Lear Ast ee he at his RG Gon ede gn Spe 2 he be Rak CR toga a es al Sa 1. 07 10. 6 120. 0 3. 68 
Libs) ae Se eS Re ReneS ce eer Caen Ne Se earn ne pA! 5. 6 22. 30 1.01 
JULY a ELS ee eR ae ee ee! SSS 33 Sao! 9. 30 . 62 
AMPUSb. 2. = seers be eget - sage os lyrobsee te ss - Peg ete 2 4E 16 | 4.9 3.85 ~ 40 
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JRA Ss LOU b ast bed h ES SR Ee es ee kee SE .10 26. 2 68. 50 4.11 
GLUT Y — eters oe of PRS dO AR OR eee on Cee . 78 31,9 793. 0 4.57 
MWienCh es... {8 OM Sew Seok ee SS? gels 2. ee eee. oo . tl 15.4 168. 0 3..00 
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INOVeItDOR eee - Me TLS Wk bee eee ee ee Sees DSi a+ seston Saal Rss ee oe Rs 
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1 Mean gage height at Dam No. 37 increased to allow for effects of local rains, as indicated by mean gage 
heights at the United States Weather Bureau gage, Madison, Ind. 
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K,H* = 131 V 


Statrons 104.-349-—— KzH*= IAG V * 
Stations 598-6iI—— KzH*= 38 V 
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RIVER AND CORRESPONDING: VALUES OF THE — 


REAERATION COEFFICIENT K, WEIGHTED BY THE SQUARE OF THE STREAM DEPTH 
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TABLE No. 6.—Relation between the reaeration coefficient (K2), the mean depth 
(H) of the stream above extreme low water, and the velocity of flow (V) 
in certain stretches of the Ohio River—Continued 





Month Ky i K, Ha Pg 
(20°. C.) (feet) second) 
Station No. 598 to station No. 611: 

Nir A014. tg id en bret Sy pre WA rt oe } 21 21, 4 96. 0 2. 94 
Ei Gp ce, Sell ainsi et oe eine faa ea eile L a i a ON ria aL a 36 6.9 Vea 85 
SUP = so 3s oo os ho ees ark corey Pa rarest 1 LT 47 6.0 17.0 68 
BRBEUELG Goto ees Se Se Nhe eS ee ee Ma ean . 43 5.1 ie? 56 
Sa RS a eens OE eee ne, DO et Senos . 87 5.5 26. 3 63 
COST OSE) Fes cee on SER ORS ys Sond Lk Oe ein rs ee a ae Seen ee . 54 5.5 16. 4 60 
INGMeI HOR secure eet oo A APTOS T AE LF SO jorres ee 39 
2 SCOP TC a Rae ee 2 nth. a RS eee em eee ee . 75 16. 2 197.0 2. 24 
dangary, LOPS. 22. si. 50 as en Ee AS BAU . 58 29.1 490. 0 3. 73 
PRS THIRRA TAY oe a ee he ae a 2. 44 36. 1 | 3, 180.0 4.11 
1 Ging 1, ee Res ee ete Meee Oe ewe ie Meee, Oe Bea . 82 16.5 223. 0 2. 28 








Notre.—Value of (H) estimated as follows for the various river stretches: 
Stations 3 to 11, from gage heights at Dam No. 3, upper gage. 

Stations 11 to 19, from gage heights at Dam No. 4, upper gage, and Dam No. 6, lower gage. 

Stations 23 to 65, from gage heights at Dam No. 8, lower gage. 

Stations 65 to 77, from gage heights at United States Weather Bureau gage, Wheeling, W. Va. 

Stations 77 to 88, from gage heights at Dam No. 17, middle gage. 

Stations 104 to 349, from average of gage heights at reference gages from Dam No. 14 to Portsmouth, in 


- glusive. 


ape 349 to 461, from average of gage heights at Portsmouth, Ohio, Maysville, Ky., and Dam No. 35, 
middle gage. 


Stations 475 to 482, from weighted average of gage heights at United States Weather Bureau gage, 
Cincinnati, and at Dam No. 37, lower gage. 

Stations 482 to 488, from gage ‘heights at Dam No. 37, lower gage. 

Stations 492 to 598, from gage heights at Dam No. 37, lower gage. 

Stations 598 to 611, from gage heights at Dam No. 41, lower gage. 


In Figure No. 8 are shown plots of (V) against (K,H?) for three 
stretches of the river selected at random, the data for these plots 
having been taken directly from the last two columns of Table No. 6. 
Tt will be noted that both ordinate and abscissa scales are logarithmic 
and that the plotted points follow closely a straight line, the equation 
of which is— 

log (K,H?) =n log V+log c. 
whence K,H?=cV® as above 


If the term (H?) be transposed to the denominator of the second 
member of the above equation, it will be noted that the value of 
(KK,) is inversely proportional to the square of the depth, as measured 
by the gage height (H). Bearing in mind that (K,) is a direct 
measure of the rate of reaeration, this empirical finding is strikingly 
in accordance with the theory of the depth-square function as a 
factor in reaeration, stated on page 28, wherein it was shown that, 
viewed from standpoints of both diffusion and dilution, the reaera- 
tion rate should bear an inverse relation to a quantity approximating 
the square of the depth. This theory is further borne out by the 
fact that when the discharge of the Ohio River, corresponding to 
(Q), is substituted for the square of the gage height in formula (4), 
the quantity (K,Q) bears the same kind of a power function relation 
to the observed velocity of flow as does (K,H?). 

In Table No. 7 are given the empirical values of (c) and (n) for 
the different river stretches, as derived from plots of the data of 


oO 


~ 


Table No. 6 similar to those shown in Figure No. 8. From the 
manner in which these coefficients have been obtained, wherein the 
biochemical factors controlling the oxygen demand and deficit. values 
have been taken into account, it is apparent that. (c) and (n) should 
be independent of such factors and should depend only upon the 
physical characteristics of the particular stream stretch for which 
they have been derived. The manner and extent of their dependence 
upon such characteristics will therefore be discussed. 


yn 
TABLE No. 7.—Values of (c) and (n) in relation: K.=¢7p, for various stretches 


of the Ohio River 
Notation: 
Ko=reaeration coefficient. 
V =mean velocity of flow, feet per second. 
H =mean depth of water above extreme low water. 















































Stretch Mean 
depth at 
ares (c) (a). } Formula 
FTO ihecl MO water 
station— | station— (feet) i 
a eat 
3 11 oP ae 0.57 | Ke=131V0-87 
He 
11 19 | 4.4 11,7 1,58 | Ke=11, 7V1-58 
H2 
23 65 bd: 2. | 6. 2 1.39 | Ko=6, 2V1-39 
a 
65 77 6.2 38. 0 1.00 | K2=38, 0V1-00 
| ; ae 
77 88 7.4 18. 0 1.48 | Ko=18. 0V1-48 
} H2 
104 349 8.8 1. 46 2.62 | Ko=1. 46V2-62 
349 461 9.5 0. 23 4.06 -| K2=0.23V4: = | 
{ H? 
475 482 12.8 60.0 | 1.71 | Ke=60. 0V!-71 
4 H? 
482 488 15,3 0. 23 §:40 | K2=0.23V5-40 
| Hi 
| 
492 598 15,2 18, 0 2, 35 | Kas. OV2-35 
| | i 
598 611 11,7 38. 0 2,12 Pores OV2-12 
| a 





1 Mean depth of extreme low water corresponds to mean depth at zero reading of reference gages. 


Values of the constants (¢) and (n) which, it will be noted, modify 
only the velocity of flow function, depend primarily upon physical 
stream conditions which influence turbulence. Broadly, the ex- 
ponent (n) defines the range of variation in rates of reaeration in 
a given river stretch under different flow conditions, while the coefti- 
cient (c) defines their general order of magnitude. Physical condi- 
tions which influence turbulence in the same river stretch under 
various flow conditions should therefore modify (n) to the greater 
extent, while those which determine the effects of the same flow 
conditions upon turbulence in various river stretches should exert 
the more influence upon (c). | 

The effect of a given variation in river stage upon the velocity 
of flow is an important factor modifying turbulence in a given 


(SW ie te ogee Se 
Os Nec aa a. ¥ 


a 


2. 





iN Feer 


STAGE 


River 


999—25. (Face p. 51.) 







FiG.No. 9 


MEAN VELOCITIES OF FLOW IN CERTAIN 
STRETCHES OF THE OHIO RIVER 


wee 


ie ! a 
Ve oa i" 
i ; 


CORRESPONDING TO VARIOUS VALUES OF ‘H” 
(Hig par SexTINe MEAN DEPTHS OF 


| 


STREAM ABOVE EXTREME LOW WATER, 





sae 


s 
- iG oe. ) 3 “ On 
gine ose os ee 
~ > - 
‘. . = 


as 


he ee 

ria ine 
ae COCO 
Pett | tt 
te P 


OL 


stretch under various flow conditions, and this effect should there- 
fore be closely related to the value of (n). The physical condition 
which influences the river stage-velocity relation is principally the 
shape of the channel. For example, channels with level bottoms and 
straight, steep side slopes produce an entirely different river stage- 
velocity relation than irregular channels with flat side slopes. In 
the latter instance it is common to observe conditions in which the 
cross-sectional area of flow of a rising river increases more rapidly 
than the discharge, with a resulting decrease in the actual mean 
velocity. In many instances in which the channel side slopes are 
straight and steep, as is frequently the case in the Ohio River, the 
velocity is closely proportional to river stage. 

The river stage-velocity relations in the various Ohio River 
stretches for which values of (c) and (n) have been obtained are 
illustrated by the curves shown in Figure No. 9, which show the 
variations of this relation in different stretches of the Ohio. In Table 
No. 8 are summarized the relative increases in velocity produced by 
specified increases of gage height in the different river stretches. 
The mean value for each stretch is an index of the river stage- 
velocity relation, and their relation to corresponding (n) values, 
as shown in Table No. 9 and Figure No. 10, is definite and fairly 
consistent, considering the comparative roughness of the index. The 
data from the stations 8-11 stretch are omitted, its river stage- 
velocity relation being very poorly defined because of the manipula- 
tion of dams during the greater part of the period of study. The 
greater part of the stretch from stations 475 to 482 is also under 
the influence of backwater from Dam No. 37 during pool stages of 
the river. For this reason the data from this stretch have been sepa- 
rated into two groups, one of which applies to “dam down” con- 
ditions and the other to “dam up” conditions. . | 


TABLE No. 8.—IJncrease in mean velocity of flow, feet per second, in certain Ohio 
River stretches as related to increase in gage height (feet) 














Increase in velocity of flow per 
5-foot increase in gage height 





T| 


tations, Ohio River stretch peers eo ee 
eights, | heights, 

5 feet to | 10 feet to| Mean 
10 feet 20 feet 











ILO Ae BA oe SE Se eo Bg plein ee ee Sie Bees PERE Une Bin cp ee Bod Wee ener 37h 
TESA aa Gee oar hg ete Phe Maat. "chr Re eI,” aie ange ire la ERM Sie ra HOI ea Bre 
Of ec hig  oRPRL — 2 a ret alin ete che oy gh Ge Ue ee eS en, eit bu ooh a aes Cad sere es eee f 6. 4 
TAS hea 20 Sytner ca a De Da ay ge, all pera Meee eat att Gro ee ese mee Oso 
OSS SA ee Sept bea We 0k PR Aa a ee Ae MEE be Bek tal ey al Cee compa daSs |e Gaede coe 1.8 
Lae Rah ie ie I Boo Ah ls les ire AM, Mia ta. la etait ele a aN ape gage hadaine eae) Noha 1.4 aly; 18 
BD see BGG AIRY TNO) — ore ge SR. ieee ge area Ase ah eaes pera e i fect he Te peeled ert ce ee? @22 
Mi Croeeae OL PEI AGH WELD) eet ee eee ee on ao RL a at etic ons ae 2.0 2. 0 
io erie Ain geen eles tpg) Eee CA a par a pS Byes a pep i 15 7 1.6 
7 Ory oe fe SOA Te See, ree geting ag Baht A Be ie a ORS A ah 1.8 2.0 1.9 
Pein alee Pry oni ae LN ee Ce a pe a AY et eh res eter a 2.6 2.0 2.3 
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TasLe No, 9.—Relation between mean increase’ in velocity of flow for given 
increase of gage height in certain Ohio River stretches and corresponding 
values of (n)* in reaeration formula (K2H’=cV"), tabulated according to 
magnitude of values of (n) é 








Ohio river stretch Mean 
velocity 
f increase 
yay of bh 5- 
atts n oot in- 
— Stations crease of 
P river 
stage 
8-a 3475-482 | 4 Negative 7.2 
4 65- 77 1.00 6.4 
3 23- 65 1.39 3:2 
5 77- 88 1. 48 6.3 
2 11- 19 1. 58 = oa § 
11 598-611 1.98 2.3 
8-b 5 475-482 2. 00 2.0 
10 492-598 2. 35 1.9 
6 104-349 2:62 1.8 
7 349-461 4. 06 as 
9 482-488 5. 40 1.6 





1 Derived from data in Table No. 8. 
2 Derived from data in Table No. 7. 
3 Dam up. 


* Probably due to influence of alge growths in October and November. Allowing for this influence, 
value of (n) probably positive and very low. 
5 Dam down. f 


As regards the coefficient (c), it has already been noted that its 
value is defined largely by physical characteristics producing differ- 
ent degrees of turbulence in various river stretches under similar 
flow conditions. The slope and irregularity of the channel are 
probably the most prominent of such characteristics. The relation 
of these two factors to values of (c) in various stretches of the Ohio 
is therefore a matter of importance. 

The numerical expression of the average slope of a river channel 
or of the water surface is simple, a convenient unit being in terms 
of feet vertical drop per mile of horizontal length. Measurement 
of the average slope of an irregular channel bottom is subject, 
however, to considerable error, since the elevation of the bottom at 
terminal points of a given stretch may not be representative, but 
may be determined by some local irregularity of profile rather than 
by the general slope of the channel. For this reason, the slope of 
the water surface at extreme low river stage is probably the more 
reliable index of the general slope of the bottom, allowing for what- 
ever change occurs in the general depth of the water. 

The irregularity of the channel, on the other hand, is not subject 
to the definiteness of expression that may be employed for the slope, 
since it is in itself an abstract, relative characteristic. Certain 
arbitrary methods of expression may be resorted to, however, which 
will permit comparison of the approximate irregularities of various 
stretches of river channel. Such methods are subject to much 


ays) 


variation, and the results derived through any one of them are to be 
regarded as merely rough indices of actual conditions. As regards 
its influence upon the turbulence of a stream, the unevenness of the 
channel is defined not only by the frequency of slope irregularities 
but also by their relative sizes and depths from the water surface. 
Obviously the expression of all of these factors, constantly varying 
in natural streams, is futile. A workable compromise which may be 
conveniently adopted consists of noting the number, per unit of 
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RELATICN BETWEEN VALUES OF “n” IN REAERATION 
FORMULA (K,H*=cV") AND INFLUENCE OF INCREASING 
RIVER STAGE UPON VELOCITY OF FLOW 


Based upon data tn Table 20 
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length, of slope irregularities, greater in size than a specified mini- 
mum, and assuming that in fairly large groups their average size 
and depth from the surface will be similar. For example, from a 


longitudinal profile of the channel] bottom the number of changes © 


in slope per mile may be counted, each one of which produces a 
change of elevation exceeding 1 foot. For convenience, this num- 
ber may and will hereafter be designated as the “irregularity fac- 
tor.” In Table No. 10, values of this factor are given for the 
various stretches of the river enumerated in preceding tables. 
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TaBLe No. 10.—Values of “ irregularity ” factor for given Ohio River stretches, 
consisting of number per river mile of changes in longitudinal slope of chan- 
nel bottom resulting in over one foot change iv elevation 








Ohio river stretch Changes 
in direc- 
fie as 

channe 

: (ode) Total bottom 
Serial Station : per mile 
No. ’ (irregu- 

larity 

factor) 

1 3-11 eA 33 4.6 

2 11-19 G3 19 | 2.6 

5 23-65 42. 2 123 2.5 

4 65-77 11.9 39 2.8 

5 77-88 | 10.8 44 4.1 

6 104-349 145. 7 375 2.6 

uf 349-461 | 127.4 286 2.2 

8 475-482 (aa. 36 4.9 

9 482-488 5.9 14 2.4 

10 492-598 105. 4 328 3. 

11 598-611 13:2 38 2.9 
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In Table No. 11 these same river stretches are divided into twe 
groups according to values of the “irregularity factor” and ar- 
ranged in each group in the order of decreasing low water slopes, 
the value of (c) being given also in each case. The first group, 
comprising stretches having factors between the limits 3 and 5, may 
be regarded as relatively uneven, while the second and larger group, 
with factors between 2 and 3, may be classified as relatively smooth. 
Both groups present, however, the fairly smooth type of channel 
found in the Ohio River and similar streams. The irregularity 
factor is taken account of in this manner, and attention may then 
be directed to a study of the relations between slope factors and 
valuey of (c) in each group. This relationship is shown in the 
table and is illustrated graphically in Figure No. 11, in which the 
characteristic curve for each group is drawn through the plotted 
values. Each relation is shown to be definite and the curves them- 
selves are consistent with-each other; for example, the curve of the 
higher irregularity group (3.0 to 5.0) gives higher values of (ce) 
with equal slopes. The. method of dividing the results into groups 
defined by limiting values of the irregularity factor is employed 
for the reason that the factors themselves are but approximate 
indices of channel unevenness and hence individual values are not 
in themselves highly significant. In streams of the type of the 
Ohio, lack of smoothness of the channel is probably not as im- 
portant a factor in reaeration as is watercourses having rough, 
irregular channels. In some of the latter type, the frequency and 
sharpness of horizontal bends is also probably a factor modifying 
the influence of channel roughness upon the turbulence of the stream. 
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RELATION BETWEEN LOW WATER SLOPE..AND VALUES OF “c” IN 
REAERATION FORMULA (K,H*=cV") ACCORDING TO ROUGHNESS 
OF CHANNEL 
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longitudinal direction of channel bottom 
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TABLE No. 11.—Relations between value of (c) in reaeration formula (K:H?= 
eV”) and low-water slope in feet per mile, for two general types of stretches 
of Ohio River, grouped according to relative irregularities of channel 


Relatively uneven ‘‘ irregularity factors” 1 





ranging from 3.0 to 5.0 
Ohio River Low- 
stretch Value | Water | “Irregu- 
a ae oF Paw Otte slope larity’”’ 
Serial Stretch (feet per | factor 2 
No. mile) 
i 3- 11 131.0 2.14 4.6 
8 | 3 475-482 50. 0 1. 69 4.9 
5 77- 88 18.0 . 50 4.1 
10 492-598 18. 0 . 21 3.1 











Relatively smooth, “‘irregularity factors”’ 
ranging from 2.0 to 3.0 





41. degetiz7 38. 0 1. 68 2.8 
11 | 508-611 38. 0 1. 67 2.9 
2 1l- 19 12.0 1. 04 2.6 
3 23- 65 6.0 65 2.5 
6 | 104-349 1.5 88 2.6 
7 | 349-461 2 36 2.2 
9 | 482-488 2 bat 2. 4 








1 “Trregularity’’ factor=number, per river mile, of changes in longitudinal direction of channel bottom- 
(See Table No. 10.) 

2 See p. 54. 

3 For open-channel conditions (Dam 37 lowered). 


From the relationship curves derived in Figures Nos. 10 and 11 it 
appears that whatever may be the weight of other factors in modi- 
fying the control exerted by stream depth and velocity of flow over 
rates of reaeration in the Ohio River (and there are doubtless other 
factors involved) the influences of channel slope and irregularity 
and of variations in river stage as related to changes in velocity are 
manifest and well defined. This is true, as will be shown later, to 
the extent that, given these three factors in the terms stated, values 
may be assigned to (c) and (n) in the formula: 

K,=" 
sufficiently close to give a very fair approximation of the true rate of 
reaeration of the Ohio under any given condition of flow. By de- 
riving empirical formule for the curves above noted, this operation 
may be done more directly. 

The relationship curve for values of (n) in Figure No. 10 is a 
hyperbola having the general form: (x-a) (y-b)=1. For this curve, 
(a)=1 and (b)=1.17, approximately. Taking the exponent (n) as 
(x) in the general equation, the following empirical formula is 
readily derived: 


1 
het bay 
or, 
HEY AAT 
Bngraon. 19 


in which (n) represents the abscissae and (y) the ordinates of the 
curve as shown in the chart. 


06 z 
e 


The two (c) curves in Figure No. 11 are more difficult to formu- 
late and unfortunately are not well fitted by empirical formule 
having exactly the same type of equation. Curve A, for smoother 
channels, has the more complex formula’ which is that of an expo- 
nential function with an additive constant, having the general equa- 
tion: y=a(10™%)+c. For the curve in question, a=0.39, b—1.16, 
and c—17, approximately. Denoting (x) by the coefficient (c) and 
(vy) by the slope (S), we then have: 


¢=0.39 (10"28) +17 = 


The formula of curve B is less complicated, being a simple power 
function of the general form: y=ax™. For this curve, a=1.1 and 
m=2.3, whence (using the same notation as above) : 


e= Gif & (6) 
Having derived formulas (4), (5), and (6), the equivalents of 
(c) and (n) may be substituted in the equation: 
ou 
1 EE 
giving the following empirical formulas for (K,) directly in terms 
of its determining factors: 
1. For relatively smooth channel (irregularity factor=2.0 to 3.0), 


K 








y—0.17 

x 10.39(10-S)417] VY (7) 
2 H?2 | 
2. For relatively irregular channel (irregularity factor=3.1 to 
5.0), ! : 
y—0.17 

1.1S237v70" 8 
K.-S (8) 


in which— 
K,=reaeration coefficient. 
V =velocity of flow in feet per second. 
H =mean river depth, in feet, above extreme low water (ordi- 


narily given by means readings of gages in or near the river stretch 
considered ). 


S=low-water slope of the channel, in feet per mile. 

y=mean relative increase in velocity of flow per 5 feet increase 
in gage height. 

The formule (7) and (8) have been derived principally in order 
to provide a concise mathematical statement of the empirical rela- 
tions modifying the reaeration coefficient (K,). While probably 
less accurate for purposes of computation than are the curves them- 
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VALUES OF REAERATION COEFFICIENT, K,, 
IN CERTAIN STRETCHES OF OHIO RIVER 
CORRESPONDING TO VARIOUS ACTUAL MEAN 
VELOCITIES OF FLOW IN FEET PER SECOND 


Based upon data in Table 23 


derived from application of general equation 
K,H*=cV" with values of “c “and “n” for each 
stretch as given in Table No. and 

values of H derived from curves in Fig.No. | 
Each value of H represents the mean a 
depth of the stream above extreme low water 
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selves, they are sufficiently well correlated with the observed data 
for all practical purposes of estimate. In applying the formulas to 


the Ohio River it is evidently preferable, where possible, to choose 
for a given river stretch the specific values of (c) and (n) which 


have been derived for it. Since, as shown in Figure No. 9, the 
gage height (H) and the velocity (V) are definitely related to 
each other, a given reading of the former in a particular section of 
the river fixes the value of the latter, permitting the direct calcula- 
tion of the reaeration coefficient (K,) for a given river-stage condi- 
tion by means of the formula: 

Tn 
K= 9 

Employing this general method, except that depths (H) were 
derived trom assumed velocities of flow (V), values of the reaera- 
tion coefficient (K,) were computed for a number of stretches of 
the Ohio River under different velocity-of-flow conditions, with re- 
sults as given in Table No. 12 and shown graphically in Figure No. 
12. In Table No. 18 are tabulated values of (H) in the several 
river stretches corresponding to the different velocities of flow as- 
sumed in obtaining Table No. 12, these values having been derived 
from the relation curves of Figure No. ¥. The purpose of drawing 
the curves of Figure No. 12 has been to show the relations existing 
between variations in flow conditions, as indicated by velocities of 
flow, and induced rates of reaeration, as measured by values of the 
reaeration coefficient (K,). 


TABLE No. 12.—Calculated values of reaeration coefficient in specified stretches 
of the Ohio River under different velocity-of-flow conditions, assuming 
Stream-depth factors H corresponding to velocities as defined by curves of 
Figure No. 9 


Value of Ke at velocities (feet per second) 
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Taste No. 13.—Mean depths (H) of Ohio River above extreme low water in 
certain stretches corresponding to given velocities of jlow 


[Derived from Figure No. 9]. 
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1 Pool stage. Dam up. 


One characteristic of these curves deserving comment is the tend- 
ency of (K,) to reach a maximum with increasing velocities of flow 
up to a certain point, beyond which lower values are obtained. The 


velocity of flow coinciding with this maximum (K,) point varies in’ 


different river stretches studied, but in a majority of them it is 
reached at somewhere between 2 and 4 feet per second. It thus 
appears that at a certain critical velocity, peculiar to a given stream 
type, the induced rate of reaeration reaches its maximum, but that 
when the velocity is greater or less than this critical value, reaera- 
tion proceeds more slowly. 

The reasons for this apparently paradoxical situation are to be 
found in the relation existing in a given river stretch between the 
two terms (cV") and (H?), in the second member of the reaeration 
formula: 

_ ev 
ac H?2 

In this formula it will be noted that if (V) and (H) increase at 
the same rate, the value of (K,) increases or diminishes according to 
whether the exponent of (V), which is (n), is greater or less than 
the exponent of (H), which is 2.0. If the exponent of (V), however, 
is less than 2.0, (K,) increases if (V) increases at a sufficiently 
greater rate than (H). The converse is likewise true; that is, (K,) 
may actually decrease in value with an exponent of (V) greater 
than 2.0 if (H) increases at a sufficiently greater rate than (V). 

Under natural flow conditions it is found (see fig. No. 12) that in a 
few stretches, as between Stations 23 and 65, the velocity (V) 
increases at about the same rate as the depth (H) above a certain 
low depth, but in a majority of the stretches (V) increases in rela- 
tion to (H) at greatly decreasing rates above a certain point; hence 


K 
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with increasing river stages the effect of depth upon (K,) becomes 
nereasingly manifest above this point, while that of velocity be- 
comes less so. In other words, at low depths the velocity effect pre- 
lominates, but at greater depths the depth factor predominates. 
The velocity-depth relations themselves are influenced almost en- 
tirely by the shape of the channel cross section, as has already been 
pointed out. 

Another characteristic of the reaeration-velocity curves is found 
in their reaching of a definite minimum point at some low velocity. 
Apparently below this point the effect of decreasing depth predomi- 
nates over that of diminishing velocity. However, these extremely 
low portions of the curves are probably not as reliable as the portions — 
ubove this minimum point; in fact, it is likely that at the minimum 
point the effect of velocity becomes almost negligible and reaeration 
tends to follow the laws of diffusion in quiescent waters. 


APPLICATIONS OF DATA 


It has been shown in the foregoing text that the oxygen self-puri- 
Heation of the Ohio River is a measurable phenomenon, governed by 
definite laws and proceeding according to certain fundamental phys- 
ical and biochemical reactions. Because of the fundamental char- 
acter of these reactions and laws, it is fairly evident that the prin- 
ciples underlying the phenomenon as a whole are applicable to virtu- 
ally all polluted streams. In considering the applications of the 
data which may be made, it is therefore important to recognize the 
relation which they may bear not only to the problems of the Ohio 
River but also to those of other streams in general. 

As regards the Ohio River, no widespread and detailed applica- 
tions of the data were justified under conditions of pollution of this 
stream, as existing at the time of these observations, for the reason 
that it was not, as a whole, seriously overburdened from a stand- 
point of its oxygen reserve. This is amply shown by the dissolved 
oxygen results which have been cited. In only two short sections of 
the river, below Pittsburgh and Cincinnati, respectively, did condi- 
tions at any season of the year approach serious oxygen depletion. 
With the possible exception of Pittsburgh, where the situation is 
complicated by a chemical withdrawal of oxygen from the river as 
2 step in the oxidation and precipitation of soluble iron salts, con- 
ditions nowhere approached an acute stage. For these reasons any 
applications which may be made of the data must deal largely with 
hypothetical cases. In this connection, however, it is instructive to 
show, by applying the principles developed from the study, approxi- 
mately what dissolved oxygen conditions may be expected in the 
river stretches below the two cities under various conditions of 
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pollution and to compare them in this respect under similar pollu- 
tion conditions. 1a oe | 
The most critical seasonal condition with respect to deoxygenation _ 
of the river consists of dry-weather flow at summer temperatures. | 
Referring to Table No. 3, the average reaeration capacities of the 
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river below Pittsburgh and Cincinnati, respectively, during the July- 
October period in 1914, are shown by the following mean values 
of the reaeration coefficient, K,: 


Below Pittsburgh, K,=0.25 (approximately). 
Below Cincinnati, K,—0.45 (approximately). 


OL 


At the mean temperatures prevailing in the river during this 
period the rate of deoxygenation below both cities is given by a 
value of the coefficient K,, equal to 0.12. 

Referring to the general formula (1) on page 18, and assuming’ 
the above constants, values of the oxygen saturation deficit (D) after 
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3 
- various times of flow (t) may be readily computed for different as- 
sumed values of initial oxygen deficit (D,) and initial oxygen demand 
(L,). Assuming these values and computing D, results are obtained 
as shown graphically in Figures Nos. 13 and 14, the former chart 
illustrating conditions below Pittsburgh and the latter parallel con- 
ditions below Cincinnati. 
999—25——5 
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In constructing the two charts it has been assumed that all of the 
pollution enters the river at a single point and is immediately mixed 
thoroughly with the stream water, producing the initial oxygen 
demand and oxygen-saturation values employed. This assumption 
is, of course, not in accordance with the actual circumstances, since 
the wastes of both cities are discharged into the river along frontages 
of several miles. In order to obtain proper values (I.) and (D,), 
it would be necessary to make a detailed study of the amount and. 
composition of the wastes discharged at each sewer outfall and also 
of the times of flow of the river between the various outfalls. 

In spite of the crudeness of the assumptions upon which they are, 
based, the two charts bring out several interesting points concern- 
ing the reaeration capacity of the river as related to its permissible 
pollution. Thus, it is indicated that no dangerously low depletion 
of the dissolved oxygen supply of the river below the two cities 
may be expected until the amounts of pollution added by them are 
sufficient to produce an initial oxygen demand in the river of 25 parts 
per million or thereabouts. Assuming the oxygen demand of nor- 
mal city sewage to be 112 grams per capita daily, this would rep- 
resent a critical dilution of 11 parts of river water to 1 part of 
sewage. The present low-water dilution at both cities 1s much 
higher than this figure. . : 

Under actual circumstances it is, of course, certain that locally 
offensive conditions may be expected to prevail in the river at both 
cities long before a general exhaustion of its reserve oxygen supply 
‘is threatened at points below them. With sluggish conditions of 
flow the rate of mixture of sewage with the river is extremely slow, 
and more or less localized zones are formed around the sewer out- 
lets within which sludge deposition and deoxygenation of the super- 
natant water proceeds far more actively than at points in the 
stream farther removed. Such zones are very likely to become foci 
of local offense, while the stream as a whole still maintains a con- 
‘siderable amount of reserve oxygen. For this reason a detailed 
study of local conditions is essential to the fixing of safe minimum 
standards of residual oxygen supply below any given city. 

As regards the application of the various constants to calculations 
of oxygen changes in different stretehes of the Ohio River, it is of 
interest to note the approximate precision obtainable therefrom in 
so far as it is indicated" by the deviation of calculated from observed 
oxygen values at the lower terminals of those river stretches for 
which data were obtained in connection with the present study. In 
order to make such a comparison values of the velocity-gage height 
factor (y) and the slope (S) in formulas (4) and (5), respectively, 
were taken for each river stretch from Tables Nos. 9 and 11, and 
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from them values of (c) and (n) were derived, using the curves of 
Figures Nos. 10 and 11. From the (c) and (n) values thus obtained 
reaeration coefficients (K,) at the stream temperature were computed 
for several river stretches and various months, following the detailed 
method outlined in Appendix C. Employing these calculated (K,) 
values, the oxygen content of the river at the lower terminal of a 
particular stretch was computed from the observed oxygen deficit 
at the upper station, employing formula (1) for this purpose. The 
calculated figure was then compared with the corresponding dis- 
solved oxygen content as actually observed at the lower station, with 
results as indicated in Table No. 14. 


Taste No. 14.—Comparison of the observed dissolved oxygen content of the 
Ohio River at various sampling stations with corresponding values calculated . 
by the use of formulae (1), (7), and (8), basing the calculations in each 
case on the observed oxygen content at the next station upstream 
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Standard deviation =a/ ee per cent. 


Probable error =0.674494.5=3.0 per cent. 
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In order to show the distribution and range in magnitude of de- 
viations of the calculated from observed values, the results in this 
table have been arranged in descending order of magnitude of the 
squared percentage deviations. It will be noted that the positive and 
negative deviations are fairly well balanced, considering the rela- 
tively small number of items. If these deviations be taken as a 
measure of the error of calculation, and if each one of the two 
groups, positive and negative, be classified according to the frequency 
of errors greater than varying amounts, a curve is obtained: which 
coincides closely with the so-called “normal” curve of error. This 
being the case, it is permissible to determine the probable error by 
the conventional method, which consists of squaring the individual 
errors (v) and computing the standard deviation and probable error 
by the following relations: 
ics 

n 





Standard deviation = a 


Probable error = 0.67449.| = 


where (2 v’) gum of squares of deviations and (n) = number of items. 


Applying this method, the standard deviation and probable error 
of the calculated oxygen values in Table No. 14 are found to be 
4.5 and 8 per cent, respectively, as referred to the observed results. 
Taking into account the errors involved in the laboratory and physi- 
cal determinations of the various factors in the formula, this com- 
posite error is remarkably low. The result obtained tends strongly 
to confirm the theory underlying the formulas applied, and to in- 
dicate that the errors involved in their application are, to a large 
extent, compensating ones. It should be pointed out, however, that 
2 calculation such as the foregoing one is by no means a test of the 
general applicability of the formule and constants derived from 
this study, for reasons that are obvious. It constitutes, on the other 
hand, an excellent check upon the formule, as far as their applica- 
tion to the Ohio River is concerned. 3 

As regards the more general application of the Ohio River con- 
stants to other streams, it is of interest to note the results of a test 
of this kind made in connection with an analysis of results obtained 
from a sanitary survey of the Thames River basin, in Connecticut, 
conducted during the summer of 1915. Among the data from this 
survey were gathered certain basic figures relative to the amounts 
of domestic sewage and idustrial wastes discharged into the Thames 
at various important points, the summer discharge of the stream 
and its tributaries, and the slope and general characteristics of its 
channel. From these data, values of the reaeration coefficient (K,) 
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were derived and, proceeding from point to point in the stream, its 
theoretical dissolved oxygen content was computed, allowing for 
the inflowing of tributaries and the introduction of fresh pollution 
at the proper intervals. The oxygen values thus calculated were 
then compared with those actually observed by the State Board of 
Health, with results as given in Table No. 15. 


TABLE No. 15.—Comparison of calculated with observed dissolved oxygen results 
at specified points in Thames River Basin, Conn., calculated values being 
based. on application of Formula (1) 
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In interpreting these results it should be pointed out that a num- 
ber of the factors entering into the calculations were virtually in- 
determinate, and broad assumptions in respect to them were neces- 
sary in order to complete the numerical computations. While there 
is considerable deviation between individual computed and observed 
values, there is unquestionably a marked correlation between them, 
as is evidenced by the fact that their coefficient of correlation? is 
0.81 and its ratio to its probable error is 11.6. It will also be noted 
that if a certain fixed standard of minimum permissible dissolved 
oxygen content, say 3.0 or 4.0 parts per million, were to be employed 
as a basis for judgment regarding the condition of the stream at 
the several stations, the calculated results would lead to about the 
same conclusions as would the observed values, as regards the zones 
in which the oxygen status of the river were satisfactory or 
otherwise. 

In general, the application to other streams of the constants de- 
rived for the Ohio River would require that particular care be 
observed in the selection of the different factors concerned ‘with 
stream and channel characteristics, which may be expected to vary 





*The Galton-Pearson coefficient. (See p. 41.) 
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rather widely from those prevailing in the Ohio. In many cases, 
us in the example above given, certain of these characteristics can 
not be properly determined in definitely measurable terms without 
detailed local study. In this phase of the matter merely a begin- 
ning has been made: in the Ohio River studies, and further work 
dealing with a large variety of stream types will be necessary be- 
fore the entire working range of reaeration coefficients can be 
established. RS 


SUMMARY AND CONCLUSIONS © 


The studies described in this paper have had as their chief aim 
the formulation of a reasonably complete set of working principles 
from which further studies of questions involving the organic sta- 
bility of polluted streams may be made with more assurance than has 
hitherto been possible. In working out such a statement it has been 
necessary, in certain parts of the discussion, to reverse the usual pro- 
cedure of developing a Jaw from empirical findings, to the extent of 
evolving what may appear to be a rather elaborate theory in order to 
lay the groundwork for interpreting the findings themselves. This 
has-been especially true of the portion of our discussion dealing with 
reneration, wherein the establishment of a basis for numerically meas- 
uring the rate factor has necessitated the formulation of a reaeration 
coefficient from wholly theoretical considerations. The basic theory 
underlying practically the whole discussion, which is summarized in 
formula (1), is in reality, however, an extremely simple one when re- 
solved into its component parts, and it is our belief that the experi- 
mental evidence on which it rests is fundamentally sound. 

The empirical relationships established with respect to reaera- 
tion rates in the Ohio River must be interpreted, for the present 
at least, as applicable only to this stream or to rivers of a closely 
similar type as regards their physiography. Further studies of 
these relations will doubtless reveal wide variations in the specific 
values of constants derived for streams of different types, and will 
probably also disclose that formulas expressing them do not always 
follow the same mathematical laws as found in the present case. 
These matters can be cleared up only by additional observations 
on streams of widely different characteristics and by the checking of 
formulas against results of observations. 

Certain basic units of measurement employed in connection with 
this study, especially those measures which apply to factors cor- 
related with stream turbulence, need far more satisfactory defi- 
nition than has been practicable in the present case. This statement 
applies in particular to the rather crude indices which have been 
employed for expressing channel irregularity and the velocity-depth 
relation. For the former, the use of the Kutter (n) or the Chezy 
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(c), and for the latter, the constants of the velocity-depth relation 
curve, might be advantageous from a standpoint of their more gen- 
eral applicability. Here again are questions for further study. 

As regards future studies of the subject, that which is probably 
most needed at present is a further experimental investigation of a 
number of fundamental points bearing on both the biochemical oxi- 
dation and the physical reaeration of streams, employing improved 
laboratory methods which have been developed since the data for the 
present study were collected. Studies of this kind should compre- 
hend the following main subjects: 

(1) Coefficients of oxidation (corresponding to K,) for various 
types of city sewage and industrial wastes. 

(2) The temperature relations of both oxidation and reaeration 
coeflicients. 

(3) Rates of reaeration of bodies of water sub jected to various con- 
_ trolled degrees of turbulence. 

(4) The influence of various physical factors, operating in natural. 
streams, upon reaeration rates, such as, for example, the action of 
dams, bends, channel obstructions of different kinds, and irregulari- 
ties in the channel bottom. The modifying influence of variations in 
stream depth upon these phenomena should also be studied under 
controlled conditions. 

Studies of this kind can be made to the best. advantage by a com- 
bination of facilities afforded by the laboratory, the experimental 
tank and fiume, and a small polluted natural stream, so situated as 
to offer opportunities for a close measurement and, if desired, a 
certain degree of artificial control, of the factors influencing natural 
purification phenomena. Further studies along these lines are con- 
templated in connection with the work of the stream pollution labo- 
ratory of the United States Public Health Service at Cincinnati, 
Ohio, where exceptionally good facilities are available for making 
them. 

In the present text the applications of the various formule and 
constants to specific problems have been merely touched upon, at- 
tention in this respect having been confined to one or two examples, 
designed to illustrate methods and to indicate roughly the expected 
precision of calculations of the kind. Further studies of this char- 
acter need to be made for a number of actual cases where the results 
of calculation may be checked against those of observation. From 
the few studies that have been made, it would appear that a high 
degree of precision is obtainable in estimating the conditions of 
pollution of the Ohio and its main tributaries under which any 
specified minimum reserve of dissolved oxygen may be maintained 
in this stream. It is also our belief that the application of the prin- 
ciples formulated in this paper to many stream pollution problems 
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of a localized nature involving the prevention of nuisance or protec- 
tion of fish life, would greatly clarify the complex points at issue 
and assist in their practical solution. With this phase of the question 
our studies have thus far been but little concerned. 

In conclusion, it may well be emphasized that the shortcomings 
of our analysis of the subject with which this paper deals are fully 
realized. Its results are presented, however, with the: hope that, 


inadequate as they may be, they may assist toward a more rational. 


treatment of the vexing problems of nuisance prevention and 
preservation of aquatic life which the alarmingly rapid increase in 
the pollution of many important rivers of the country has entailed. 
These problems are largely economic ones, and in solving them, full 
account should be taken of the enormous asset represented by the 
reaeration capacities of flowing streams. This asset is not only tangi- 
ble but definitely measurable, and it represents a saving of millions 


of dollars to the cotintry at large if properly utilized and rightly 


distributed. As regards the most advantageous plan of its utiliza- 
tion, and more especially of its just distribution, an engineering and 
legal problem is thereby presented which may well challenge the best 
thought among, those charged with the public welfare. 





APPENDIX A 
| 
SOLUTION OF DIFFERENTIAL EQUATION GIVING RATE OF CHANGE IN DISSOLVED 
OXYGEN WITH COMBINED REAERATION AND BIOCHEMICAL OXIDATION 
Notation: 


D,=initial oxygen saturation deficit in parts per million. 

D= oxygen saturation deficit at time (t) in parts per million. 

K,=deoxygenation coefficient. 

=reaeration coefficient. 
t= time in days. 

L,= initial oxygen demand in parts per million. 

Differential equation: . 
a LK D | 
= K,L,e-mt— K,D (1) 

This equation is linear, first order (Leibnitz), of general form— 


dy = 
| ¢ Aah Py=Q 
(see Murray’s “‘ Differential Equations,” 


and Q=K,,L,e—*t. 
Let iB . 


Ther = 4+K:D= O 


oe 


pp. 26 and 27) where y=D, x=t, P= Ky, 


dt | 
Log D= — f Kedt 
D= Ce—Sk2t = Ce—Est 
DeSsk:dt = C (2) 
Differentiating (2) 
eS kK:dt(d D-+ K,Ddt) = 
eS k2dt (qd D+ K,Ddt) = 8 ins 1a e—kit dt 
DeSsk:dt= fe/k:dt K,L,e-*¥:t dt+ C 
D=e—S xt (fe/k:dt K,L,e-*t dt+C) 
=e-Ket (K,L, fex2t e—kit dt+C) 
=e—Ket (K,L, fee-knt dt+C) 


—e—Kot ete emt) + Coral 








(3) 
In order to check (3) at this point, differentiate it thus: 
dD KY La ze 
eS ee e—kEit — CK.oe— Kat 
do Rie Gk, oan kw dears 
Multiplying (3) by (IKk,). 
Danae e-Kit + CK,e—x:t (5) 
Adding (4) and (5): 
dD Bihan AK) tee 
rie - KD =—_ mas: L,e-*% 
4 ates ae (this is same as (1)). 
To determine the value of (C) from (3), let “47 =O 
* KL 
Then ae ie Ke 
and C= D.- itn 


a ee K, 
(69) 


Then (3) becomes 





K,L, KyL ; 
pe ia | 
D= oo eS Kat + D,e—Kat (6 
— Bibs (e—kit — e—Ket) +. D.e—Kat (7) 
ko — Ky oe 





APPENDIX B 


Reformulation of resultant oxygen formula (1)! for a condition such that any 
deviation of (Lp) from (Lae-*:t) is assumed to be due to pollution or dilution 
uniformly distributed along the stream between stations A and B. Comparison 
of results obtained by application of the equation thus developed, with results 
obtained by using formula (1).! 


Given: 

La =the observed oxygen demand at station A. 

Lp =the observed oxygen demand at station B. 

L’,==the oxygen demand at station B as calculated from the observed value 
of (La) at station A. 

t =time of flow, in days. 

K, =—coefficient of deoxygepation. 

K, =coefficient of reaeration. 

L, increase or decrease in oxygen demand produced up to time (t) by uni- 
formly distributed increments of inflow. 

p increase or decrease in oxygen demand of stream, per unit of time, pro- 
ducing an effect equivalent to (Ly) at time (t). 

At station B, (Li) =(Ls—L’,). 
Required: 


Integrated form of differential equation oP = K, L—K.D for the condition 


assumed, 
where D=oxygen saturation deficit at time (t) 
and L=remaining oxygen demand at time (t). 
For the assumed condition, . 
L= Lae-mit+ Ly 


or a= ik (Lae—kit+ L,) — KD (1) 


The problem first becomes one of determining the value of (L;) in terms of de- 


terminable factors. In order to derive an expression for (L;) in terms of (p), 
-the method of increments may be employed, assuming (L,) to remain constant 
through each unit increment of time (At). 
Then, where t=1, Li=p 
where wey L.=p+pe-— 
where t=3, L,=p+pe-F:+ pe—2k 








where t=n, Lt=p-+pe—*1-++ pe—8Ki_. ___- pe-@—)) ki 
From this series the sum of the terms up to any time (t) may be derived thus: 
2bO es) 

hails BASE A (2) 

—e-kit 

at Station B, BY ane a ee Ae i) 
—e-Ki 
. (Lg — L’,4) (1—e7-*1) ; 

or, p= GB (3) 











| 


————? 


fe 
Substituting the value of (L,) from (2) and (1), 
ee ae ee ee 2 |- KD (4) 


This lie is a linear one of the first order, corresponding to Leibnitz’s: 


aa oY 4 py == () where x==t 





y= 0 
P= Ke ; 
Lp 
Q=T,e-mt PGE. 
Le —Kit 
Let | Tyemnt Po Ce |=0 
dD i 
Then = rabrinal 
and De/sk:dt==¢ (5) 


Differentiating (5), we have 
efk2dt (dD-+ K,Ddt) = 
p (1— e- Kit) 
=eJ/Kidt, Ky (Tse-mt P= a Sd 


t 
DeSxxdt= feskdt, K, (Dre PES Gd pte = ‘)at-+C' 


D=e-Smidt fesmdt, K, (Tyemty PSE 
KilLada— —e-*1) — —K,P KyP | Pp aliest 
—Kot ==KaNt NS sialon int |) 
(K,—K)(—e-*) KG —e= 4) Te Sige sae 
To determine the value of (C’), let (t)=0in (6). Then (6) becomes 
~~ 
Ky, Ls (1—e- Ki) — Kip 


a (KK, — Ky) (1—e-®:) pee as (i ea) eae (7) 


)at+ Cle sxaat 


=e 








Now from (3), 


{Lp ~ b's) (1—e7*) 


p= i e— rt 


Substituting the equivalent of (p) into (7) and canceling out the terms 
(l=e,), we have at Station B, 











ee ee 2 oer 
ts =e aire < = : 1 ere nt (e—Ket — e—Kat) -}- eS i as 
(1—e—kst) + Dae—Est (8) 
KL K,(Lg— LV’ 
= Gre, (et em + Demet | RUPE RES (1 —e-mit) 
A iy be. Bane ap A) Sa aes 
GE) K) a i= Erk) (e Kit —e xt) | (9) 


Referring to the original formula (1), page 18, it is noted that equation (9) 
differs from it only in respect to the rather complex expression contained in the 
_ brackets, [], which expression represents the effect of the assumption of 
uniform distribution of inflow between stations “‘A’’ and ‘‘B.”’ 
Simplifying further equation (9) by substituting for (l’,) its equivalent 
(Lae—kit), we have 
Kyl, 


“A 7h a lire tec RR al Sala 


Ki (Lp — Lye—it) cay Ki (ip— Lye st) ‘S 
+ Reena gc ees em eres) | (0 
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Equation (10) reduces to the symmetrical form: 


_ Ki Ko(1—e—mt) (Lp — Let) — Ki (Ln Lge") de®) yt (11 
ae KK (Ko— Ky (1—e-*) +Dae-=t (11) 


Comparison of results obtained in calculations of (K:.) for four stretches 
of the Ohio River, by use of formula (1), as applied in the main text and by 
use of equation (11), as developed by assuming a uniform distribution of inflow 
between stations A and B: 

In order to test the application of the original resultant oxygen formula 


(1), employed in connection with the main text against results obtained by 


assuming that any deviation of the observed oxygen demand (Ly) at Station 
B from the Station A residual (L’,) at Station B is due toa uniformly dis- 
tributed inflow between the two points, a comparison was made* of (Ke) 
values obtained in four stretches of the Ohio River, first by using the original 
formula, (1) and, second, by employing equation (11), developed in this 
appendix. 


To make the test severe, the longest two river stretches of the entire series — 
(viz, stations 104-349 and stations 492-598) were included in the four 


stretches selected for comparative study. The other two stretches selected 
(stations 475-482 and stations 77-88) are of less than average length; hence 
the two extremes of stretch length are fairly represented. The data for the 
calculations were taken from Table No. 2, presented in connection with the 
main.text. The method of calculation employed was similar to that which 
is outlined in Appendix C, modified in applying it to equation (11) only in 
respect to the terms used. (See derivation of equation (11), this appendix.) 

A summary of values of (K:) obtained by use of the two formule is given 
in the table following. 





F Value of | (K:) by | Devia- | Per cent 









































formula |equation| tion devia- 
River stretch : Month () (1) (a) tion 
from referred 
(a) (b) (b) to (b) 
1914 
Stations 104-349 (245 miles) .-.----------------- 1h (>, aa Se 0. 180 0.165 | +0. 015 9.1 
JUNO ee os . 268 . 299 —. 031 10. 4 
h Pasiy Gs. . 202 188 | +.014 7.4 
AeIStius 5-eS 212 . 190 +. 022 1L.6 
September - --- fed ep) SAb2 +. 020 13:2 
October-___.--- . 200 . 187 +, 013 6.9 
NA GAY Fo eel oS See SUS eet ee ain | oe a Sa ae Ee oe |-~---~-~--|]------=--- 9.8 
| 1914 
Stations 492-598 (106 miles) _-...-.---- June_s ke baz . 708 —.176 24.8 
Julivec eee 1. 090 1, 792 —. 702 39. 2 
Aumist 2.35.2 2s . 788 B16 +, 273 53.0 
September ---- . 761 1. 084 —. 323 29.8 
October =------ . 247 we —,. 029 10.5 
Wéan.... nos enna ch 2 oda ek endow nee =<] op Bee at eet eee FEF |e oo a een 32. 4 
1914 
Stations 470-482 (7 mules)... ae sae sk ci ape . 225 206 +.019 9.2 
JiMesse 2b See . 570 575 —. 005 9 
UL yeeee ene . 650 BOOT —.007 1.1 
AUTUStHI oO: « . 545 . 550 —. 005 .9 
September - --- . 514 . §12 +. 002 4 
Ocbobers sone a2 . 601 . 600 +. 001 a | 
1915 
IMSS co eee 802 . 805 —. 003 4 
VUNCs eae eee 1. 005 1. 056 —, 051 4. 
Meatl...u oon ne nnn we cw enna ne we nme ww nena ws lone eecncanenn--[3si--75- hs] 2--98-- Sa] ence esaene 





' 


2 The calculations incidental to this comparison were made by Sanitary Engineer H. R. 
Crohurst, of the Public Health Service. 
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Value of | (K2) by | Devia- | Per cent 





formula | equation tion devia- 
River stretch Month (1) (11) (a) tion 
from referred 
(a) 0b) (b) to (b) 
ais 529 541 012 2, 2 
Stati Ave Ss. Gatsmmiles) i “tte sess tat eee drahyetsecck 3 , -. ; 
hace ; AAI SUSt oo 2S . 245 . 236 +. 009 3.8 
September - --- . 716 . 784 —. 068 8.7 
October_.--.-- . 144 . 136 -+. 008 5.9 
Iiyifcle) Pep ee Meee keen SA EMER ELD pe (ame Gees Gee SNe aS Sige oe SPOR I || See ee roe lee rere nl eee ge 5. 2 


In interpreting the results above tabulated, it should be pointed out that, for 
purposes of application, values of (K2) are not to be considered as being sig- 
nificant beyond the second decimal place; hence deviations in the third decimal 
place are insignificant. It will be noted further that the deviations obtained in 
the present series of comparisons, when taken as a whole, are fairly well 
balanced as regards algebraic sign, 11 of them being positive and 12 negative. 
This would appear to indicate that the differences in (K:) values as obtained 
by the two methods are not due to any inherent source of error in one pro- 
cedure or the other, but are more likely due to departures, in individual cases, 
of conditions of inflow and other disturbing factors from those which may be 
- assumed in applying a particular formula. That the errors involved in either 
procedure are not serious as far as the final result is concerned is indicated by 
the close agreement of (K.) values obtained by the two formulae in three out 
of the four cases for which the comparison was made. In one of the three 
cases where the agreement was reasonably close (viz, the river stretch, stations 
104-349) the degree of concordance between the two sets of figures is surpris- 
ingly high in view of the extreme length of this river stretch and the corre-, 
spondingly long time period over which the calculation must be made. 

' In the single instance in which the agreement between the calculated values 
of (K:) was not close (viz, the stations 492-598 stretch) there is good reason 
for believing that the use of equation (11) would give results subject to a 
considerable margin of error, since a large proportion of the inflow entering 
the river stretch in question is discharged into it at a single point, about mid- 
way between the two terminal stations (through the mouth of the Kentucky 
River). An assumption of uniformly distributed inflow in a case of this kind, 
such as would be involved in the use of equation (11), would be an obvious 
source of error, tending to give higher than the true figure. On referring to 
the table it will be noted that in four of the five months for which the com- 
parison was made for the stations 492-598 stretch the values of (K:) obtained 
by use of equation (11) were higher than corresponding values obtained by 
formula (1). | 

With the one exception above noted, the differences in values of (K:) as 
obtained by the two formulae are not strikingly great, their general agreement 
being well within the expected limits of precision in calculations of this kind. 
Furthermore, the evidence at hand does not indicate the existence of any 
fundamental source of unbalanced error in the use of formula (1), such as 
would vitiate it for purposes of caleulation where conditions of disturbance in 
a given river stretch are indeterminate. For such cases, as well as the one 
most frequently encountered in practice, namely, a condition of massive pollu- 
tion of a stream at certain well-defined points, the wisdom of employing for- 
mula (1), rather than the much more complex formulae of the type of equation 
(11), appears to be borne out rather strikingly by the results obtained from 
the foregoing comparison. 


v 
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APPENDIX C 
METHOD OF DERIVATION OF REAERATION COEFFICIENTS FOR OHIO RIVER 


A. Data.—(1) Monthly mean initial dissolved oxygen content at upper. and 
lower ends of stretch. 

(2) Monthly mean final dissolved oxygen content at upper and lower ends of 
stretch, after 24 hours’ incubation in closed bottle at mo. (Other incubation 
periods=2 or 5 days.) 

(3) Monthly mean 24-hour loss of dissolved oxygen upon theghation at 20° 
C., obtained as differences of monthly mean initial and final dissolved oxygen 
content. (Other incubation periods=2 or 5 days.) 

(4) Temperature of water at each end of stretch (monthly ae 

(5) Time interval of flow for stretch, in days (monthly mean). 

(6) Mean velocity of flow for stretch in feet per second (monthly mean). _ 

B. Method of calculation—(1) Using value of 9=1.047, correct 24-hour 


losses upon incubation from 20° C. to T° C. (river temperature), assuming | 


kX, =0.1 for 20° C. 
(K,) is value in formula: Log = K,t 
L’= initial total oxygen demand. 
L =final total oxygen demand after time t. 
(2) Referring to Figure No. 6 b, compute the total oxygen demand (L) 
from the observed 1, 2, or 5 day loss of dissolved oxygen on ineubation, at 
20° C., of samples pee a Station A. If the incubation temperature is 20° C. 


PO D6 2 1F596.969 + eae & 





Ae l-day loss 2-day loss 5-day loss 


(3) Compute the total oxygen demand (L,) at Station B by a similar 
procedure. 
(4) With the value of (Ki) at the river temperature, compute the value of 
(La) thus: : 
py bat balk) 


(5) Formula for (K,) : 


Dea gi 138 eR 10—Ket) +D, - 10—Kat 

Da—dissolved oxygen saturation deficit at Station A (parts per million). 

D»= dissolved oxygen saturation deficit at Station B (parts per million). 
~ La=mean total oxygen demand (parts per million) at Station (A), based on 

mean of (Lg) and (Ly). 

K,.—deoxygenation coefficient at river temperature. 

t =time from (A) to (B), in days. 

K.=reaeration coefficient at river temperature (unknown). 

(2) Assume three or four values of (K:) and solve for (Dp) in each case. 

(3) Plot values of (K:.) and (D») obtained as above, and connect points with 
smooth curve. 

(4) Take from the curve the value of (K.) corresponding to the actual value 
of (Dp). 

(5) Substitute the value of ( K:) obtained as in (4) in the formula and com- 
pute (Dv). If (K:2) is correct, (Dv) should be the same as the observed value. 

Correction of (K:) values for temperature: Use temperature relation curve, 
-as given in Figure 7, correcting (K:) from the river temperature to 20° C, 





(i) 


Correlation of (K:) at 20° C. with mean velocities (V) and depths (H): 


n 
(1) General formula: Ka 





where (V)=mean velocity in feet per second and (H)=mean depth of the river 
above extreme low water, taken in the present case as the gauge height at a 
convenient reference gauge for the river stretch in question. 

(2) Plot corresponding values of (K:H’) and (V) on logarithmic paper, with 
(K2) as ordinates and (V) as abscissae. 

(3) Seale the slope of a line drawn through the plotted points. This slope, 
expressed as a ratio, equals the term (n) above. 

(4) Determine (c), which is represented by the intercept on the vertical axis 
of the plot. 

(5) The equation of (K:) in terms of (V) and (H) is now defined, and the 
method of correlating (c) and (n) with physical conditions has been given in 
the main text. ~ 

C. Numerical example— 

Station A=Station 3, Ohio River 

Station B=Station 11, Ohio River 

Mean temperature (T) =24.6° C. 

Time of flow (t)=1.81 days. 

(K,) at 24.6° C.=0.124. 

D,=8.36—4.36=4 parts per million. 

D,=8.53—7.21=1.32 parts per million. 

L,=2.73 parts per million. 


® 


\Month of July, 1914 (Table No. 1). 


Ly Ky 0.885 


107 ®t 9.598 





Dp = 0.89 X0.474+ 
Try Kj=-0.500 10° *"*=0.598 + 4.00 0,.124= 


K,=0.124 10,124 _ 0-385 5 gg 0:422 x 0.496 = 
K,—K,=0.376 diff. =0.474°~0.376 °°” =0.918 
This value too low. Therefore (K.) assumed too high. 
Dg=1.90 0.3812 


Try K,=0.300 107***-~0.598 X 4.00 0.286 = 
K=-0124 . 10° *"—0286 —- 90-385, 99 = 0-594+1.144= 
K,—K,=0.176 diff. =0.312 “0.176 °"" =1.738 
This value still too low. 
Dg=4.41 X0.163+ 





Try, k,—0,200 20°" — «508 + 4.00X0.435= 
K,=0.124 107**—0.435 0.385 _4 4, =0.718+1.74= 
K,—K,=0.076 diff. =0.163" 0.076. ~~ =2.46 


Plot above corresponding values of (K,) and (Dg) and take off value of (K,) 
eorresponding to Dg=1.32. Thus, K,=0.38. 


| To confirm K,=0.38. 


Try K,=0.380 107**—0.598 Dg= 1.31 X392+ 
+4.00X 0.206 = 
K,=0.124 107**=0.206 ,_0,335_1 3, =0.51+0.81= 
K,—K,=0.256 diff. =0.392 ER 6 =1.31 
This value checks that of (Kz) fro cee 
VIBRARY 


i 
Ve eT AY 


















